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We study numerically and analytically the behavior of vortex matter in artificial flow channels 
confined by pinned vortices in the channel edges (CE's). The critical current density J 3 for chan- 
nel flow is governed by the interaction with the static vortices in the CE's. Motivated by early 
experiments which showed oscillations of J 3 on changing (in)commensurability between the chan- 
nel width w and the natural vortex row spacing bo, we study structural changes associated with 
(in)commensurability and their effect on J 3 and the dynamics. The behavior depends crucially 
on the presence of disorder in the arrays in the CE's. For ordered CE's, maxima in J 3 occur at 
commensurability w — nbo (n integer), while for w 7^ nbo defects along the CE's cause a vanishing 
J 3 . For weak disorder, the sharp peaks in J 3 are reduced in height and broadened via nucleation 
and pinning of defects. The corresponding structures in the channels (for zero or weak disorder) 
are quasi-lD n row configurations, which can be adequately described by a (disordered)sine-Gordon 
model. For larger disorder, matching between the longitudinal vortex spacings in and outside the 
channel becomes irrelevant and, for w — nbo, the shear current J s levels at ~ 30% of the value J 3 
for the ideal commensurate lattice. Around 'half filling' (w/bo — n ± 1/2) the disorder leads to new 
phenomena, namely stabilization and pinning of misaligned dislocations and coexistence of n and 
n ± 1 rows in the channel. At sufficient disorder, these quasi-2D structures cause a maximum in J 3 
around mismatch, while J 3 smoothly decreases towards matching due to annealing of the misaligned 
regions. Near threshold, motion inside the channel is always plastic. We study the evolution of static 
and dynamic structures on changing w/bo, the relation between the J 3 modulations and transverse 
fluctuations in the channels and find dynamic ordering of the arrays at a velocity with a matching 
dependence similar to J 3 . We finally compare our numerical findings at strong disorder with recent 
mode-locking experiments, and find good qualitative agreement. 

PACS numbers: 74.25.Qt flux motion; 83. 50. Ha flow in channels ;83.50.Lh slip boundary effects ; 71.45.Lr 
CDW's; 62.20.Fe Deformation and plasticity 
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I. INTRODUCTION 



FIG. 1: Sketch of the artificial channel geometry. In the gray 
areas vortices are pinned by the strong-pinning NbN layer 
27 while inside the channels pinning due to material inhomo- 
geneities is negligible. The etched channel width w et ched (of 
the order of a few row spacings &o) and the effective width w 
are indicated. 



The depinning and dynamics of the vortex lattice (VL) 
in type II superconductors is exemplary for the behavior 
of driven, periodic media in presence of a pinning poten- 
tial 0. Other examples range from sliding surfaces ex- 
hibitingstatic and dynamic friction and absorbed mono- 
layers p| to charge density waves (CDW's) Jj, Wigner 
crystals and (magnetic) bubble arrays jj- Vortex 
matter offers the advantage that the periodicity clq of the 
hexagonal lattice can be tuned by changing the magnetic 
induction B. In addition, the effect of various types of 
pinning potentials can be studied. This pinning poten- 
tial, arising from inhomogeneities in the host material, 
can be completely random, as in most natural materi- 
als, or can be arranged in periodic arrays using nano- 
fabrication techniques 0] ■ In a variety of cases corre- 
lated inhomogeneities occur naturally in a material, such 
as twin boundaries and the layered structure of the high 
T c superconductors 

Depinning of the VL in a random potential gene rall y in - 
volves regions of plastic deformations 0, ITolllll IT^.Il3| , 
i.e. coexistence of (temporarily) pinned domains with 
moving domains. For very weak pinning the typical do- 
main size can exceed the correlation length R c of the VL 
(see 03) and the weak collective pinning theory can 
be successfully used to estimate the critical current den- 
sity J c 0, 0| . However, as either the ratio of the VL 
shear modulus cqq and the elementary pinning strength 
or the number of pins per correlated volume decreases, 
plastically deformed regions start to have a noticeable ef- 
fect on J c . Recent imaging experiments |l7j have shown 
directly that the rise in J c in weak pinning materials 
near the upp er critical field B C 2, known as the peak ef- 
fect llrJ ]. originates from such, rather sudden, en- 
hancement of the defect density. This strong reduction 
of the VL correlation length is also accompanied by a 
qualitative change in the nature of depinning: for strong 
pinning, depinning proceeds through a dense network of 
quasi-static flow channels (filaments) such that the typ- 
ical width of both static and moving 'domains' has ap- 
proached the lattice spacing 

dEl Ell!!- Depinning 
transitions via a sequence of static, channel like struc- 
tures have also been observed experimentally via trans- 
port experiments |2*oj . 



In superconductors with periodic pinning arrays 
(PPA's), matching effects between the lattice and the 
PPA become imp ortant. As shown first by Daldini and 
Martinoli |2lll22l |. when the vortex spacing coincides with 
the periodicity of the potential, pronounced maxima can 
occur in J c , while at mismatch defects (discommensu- 
rations) appear which gives rise to a reduced J c . The 
last decade, many more studies of VL's in PPA's have 
appeared, both experimentally and numerically. Pro- 
nounced commensurability effects were found in films 
with 2-D periodic pinning jo, HE! f° r nux densities equal 
to (integers of) the density of dots. In these systems, vor- 
tex chains at interstitial positions of the periodic arrays 
(e.g. at the second matching field of a square pinning 
array) can exhibit quasi one-dimensional motion under 
the influence of the interaction with neighboring, pinned 
vortices (24 , as has also been observed in numerical sim- 
ulations p5| . In addition these simulations have revealed 
that, depending on the vortex interactions and the sym- 
metry or strength of the PPA, a rich variety of other 
states and dynamic transitions can occur, often leading 
to peculiar transport characteristics. 

Besides the above examples, the phenomenon of vor- 
tex channelling can also arise from the presence of grain 
boundaries in the sample. Historically, the 'shear' depin- 
ning of vortices in grain boundaries in low T c materials 
received considerable attention pfj l27j because it could 
explain the quadratic decrease of J c near B C 2 in prac- 
tically relevant poly-crystalline superconductors. More 
recently, the interesting issue of channelling of mixed 
Abrikosov- Josephson vortices in grain boundaries in high 
T c materials was addressed in detail by Gurevich (2^ • 

A system in which channel motion and its dependence 
on the structural properties of vortex matter can be stud- 
ied systematically is that of narrow, weak pinning flow 
channels in a superconducting film [22], see Fig. ^ The 
samples are fabricated by etching straight channels of 
width w e tched ^ 100 nm through the top layer of an a- 
NbGe/NbN double layer. With a magnetic field applied 
perpendicular to the film, vortices penetrate both the 
strong pinning NbN in the channel edges (CE's) and the 
remaining NbGe weak pinning channels. The strongly 
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FIG. 2: Data: critical shear force density F s = J S B, deter- 
mined using a velocity criterion v/ao ~ 1 MHz, versus ap- 
plied field for a channel sample with w e tched ~ 230 nm at 
T = 1.94 K. Drawn line: Eq.Q with A = 0.05, B c2 = 1.55 
T, X(T) = 1.13 fim and effective width w = 300 nm. 

pinned CE vortices provide confinement to the vortices 
inside the channel as well as the pinning (shear) potential 
which opposes the Lorentz force from a transport cur- 
rent J applied perpendicular to the channel. By chang- 
ing the applied field H one can tune the commensurabil- 
ity between the vortex lattice constants and the channel 
width, allowing a detailed study of the shear response and 
threshold for plastic flow as a function of the mismatch 
and the actual microstructure in the channel. 

Phcnomcnologically, plastic flow in the channel occurs 
when the force density F = JB (with B/$ the vortex 
density) exceeds 2T max /w, where r max = Acqq is the flow 
stress at the edge (the factor 2 is due to both CE's) and 
w is the effective width between the first pinned vortex 
rows, defined in Fig. ^ Thus, the critical force density 
is given by: 

F s = J S B = 2Ac 66 /w. (1) 

The parameter A describes microscopic details of the sys- 
tem: it depends on lattice orientation, (an)harmonicity 
of the shear potential, details of the vortex structure in 
the CE's and the microstructure of the array inside the 
channel. Critical current measurements as function of ap- 
plied field reflected this change in microstructure through 
oscillations of F s , shown in Fig. for a channel with 
Wetched ~ 230 nm. Note that in such a narrow channel, 
the pinning strength due to intrinsic disorder in the a- 
NbGe is at most 10% of F s (except for B < 50 mT) and 
does not affect the oscillations. Since the natural row 
spacing of the VL is b = \/3a /2, with a% = 2$ /v / 3~B, 
and in our geometry B ~ u H one can check that the 
periodicity of the oscillations corresponds to transitions 
from w = nbo to w = (n ± l)bo with n integer, i.e. the 
principal lattice vector a*o is oriented along the channel 
(Fig. The envelope curve represents Eq.Q with 

Brandt's expression for the VL shear modulus |30l | : 

c 66 = ®° Bc2 M l - b) 2 (l - 0.58b + 0.29b 2 ), (2) 



(b = B/B C 2 is the reduced field and A is the penetra- 
tion depth) and a value A — 0.05. This value for A 
is close to the value J (u 2 ) /ao = 0.047 for the relative 
displacements at the crossover from elastic to plastic de- 
formations as obtained from measurements on the peak 
effect 0, ^3 ■ This led to a qualitative interpretation of 
the reduction of F s at minima as being due to defects in 
the channel, which develop at incommensurability. How- 
ever, recent developments 0, 152, H3, 13 have shown 
that (strong) structural disorder may be present in the 
CE arrays, in which case the interpretation can drasti- 
cally differ. 

In this article we present numerical and analytical 
studies of the threshold force and dynamics of vortices 
in the channel system for various degrees of edge dis- 
order. In an early paper [35j we studied the commen- 
surability effects in the idealized case with periodic ar- 
rays in the CE's. In this situation F s at matching 
(w = nbo) is equal to the ideal lattice strength 2A°cqq/w 
(the value A = 1/ (ttvo) follows from Frenkels consid- 
erations |36|). while at mismatch dislocations develop, 
leading to A ~ 0. The resulting series of delta-like peaks 
in F s versus matching parameter differed considerably 
from the experimental results, which could not be ex- 
plained by thermal fluctuations or intrinsic disorder in- 
side the channel. Therefore we investigated the effect of 
positional disorder in the CE arrays on F s near commen- 
surability (w f=a nbo) H3- I n this regime the behavior is 
dominated by the longitudinal displacements of vortices 
in the chains, i.e. quasi one-dimensional (ID), and F s is 
controlled by defects with Burgers vector along the chan- 
nel. At weak disorder, we found a clear reduction of F s 
at commensurability caused by nucleation of defects at 
threshold, while the existing defects at incommensurabil- 
ity become pinned by disorder, leading to an increase of 
F s in the mismatching case. 

The present paper first describes in detail these quasi 
ID phenomena near commensurability and/or for weak 
disorder. Using a generalized sine-Gordon model, we 
quantitatively describe how the structure and transport 
properties depend on the vortex interaction range and 
on weak disorder in the CE's. Besides the connection 
to our system, these results also provide a background 
for understanding quasi-lD vortex states and matching 
effects in artificial PPA's, including the effects of disor- 
der which these PPA's may contain due to fabrication 
uncertainties. 

The ID model shows that, above a certain disorder 
strength, spontaneously (disorder) induced defects along 
the CE's dominate over incommensurability induced de- 
fects. The commensurability peak in F s is then com- 
pletely smeared out with a value of F s at matching 
(w = nbo) saturating at ~ 30% of the ideal lattice 
strength. In the more general case of wider channels, 
the transverse degrees of freedom, especially away from 
matching (w/bo ~ n±l/2), lead to new phenomena: un- 
der the influence of disorder, the channel array may split 
up in regions with n and nil rows, involving disloca- 
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tions with Burgers vector strongly misaligned with the 
CE's. At sufficient disorder strength, such dislocations 
lead to a more effective pinning of the array then the 
'aligned' dislocations around matching. F s then exhibits 
a smooth oscillation as function of w/bo, similar to Fig. |2 
with yield strength maxima occurring around mismatch. 
This behavior resembles the classical peak effect, i.e. at 
mismatch the enhanced ability of the arrays to sample 
configuration space allows better adjustment to the ran- 
dom CE's. In the last part of the paper we show detailed 
simulations of both static and dynamical aspects of this 
behavior, including a study of reordering phenomena at 
large drive. We find an ordering velocity of the arrays 
with a channel width dependence similar to that of the 
threshold force. U sing a modified version of the dynamic 
ordering theory in |3jj it is shown that such behavior can 
be explained by a reduction of the energy for formation of 
misaligned defect pairs away from matching. The numer- 
ical results at strong disorder are also in good qualitative 
agreement with recent mode-locking experiments on the 
channel system 0, . 

The outline of the paper is as follows. In Sec[n] the 
channel geometry and the simulation procedure are dis- 
cussed. The first part of the paper deals with channels 
having hexagonal, ordered arrays in the CE's: in Sec lIIII 
we present the sine- Gordon description and numerical re- 
sults for a single ID vortex chain in an ordered channel; 
in section llVl wc show how the ID behavior extends to 
wider channels with multiple rows and ordered CE's. The 
second part of the paper deals with disordered channels: 
sectionlyldescribes the effects of weak CE disorder on the 
behavior of a ID chain, both analytically and using nu- 
merical simulations. The effects of weak disorder in wider 
channels are discussed in Sec lVIl Section TVIII describes 
the static and dynamic properties of wider channels in 
presence of strong disorder, including an analysis of the 
reordering phenomena in this situation. A comparison 
with the dynamic ordering theory, the confrontation with 
experiments and a summary of the results are presented 
in Secs lVIlTl andlrXl 

II. MODEL AND NUMERICAL PROCEDURE 

We consider straight vortices at T = in the geome- 
try as illustrated in Fig. [3| for the case of 1 row in the 
channel. The approximation T = is well justified over 
a considerable range of the experiments (see Sec lVIIl|l . 
The CE's are formed by two semi-infinite static arrays. 
The distance between the first vortex rows on both sides 
of the channel is w + bo, with w the effective channel 
width. The vortices are assumed to be fixed by colum- 
nar pins in the CE's. The principal axis of the pinned 
arrays is along the channel direction x. A relative shift 
Aa; is allowed between the arrays. In (a) the simplest 
configuration is shown, where CE vortices form a perfect 
triangular lattice. For Aa; = 0, their coordinates are: 

r„, m = ([n + frac(m/2)]a , m[b + (w - 6 )/2|m|]), (3) 
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FIG. 3: Channel geometry with pinned vortices in the gray 
areas. The specific case of 10 ~ bo, i.e. with 1 row in the 
channel is illustrated, (a) Ordered situation: the equilibrium 
positions r n , m of pinned vortices in the CE's are denoted by 
(o). The effective width w and relative shift Aa; are indicated, 
(b) Disordered case. The randomized vortex positions are 
denoted by (•). 

for ?7i =/= and frac(m/2) denotes the remainder of m/2. 

Disorder is incorporated in the model by adding ran- 
dom shifts d to the coordinates of the ordered arrays: 

I^n,m — I"n,m ~t~ d n m , (4) 

as shown in Fig. EHb). The amplitudes of the ran- 
dom shifts are characterized by disorder parameters A^, 
and A y as follows: transverse relative displacements 
d v /ao are chosen independently from a box distribution 
[—A y , Ay]. The longitudinal shifts d^ m are chosen such 
that the strain (d n+ i — d n )/ao along the channel is uni- 
formly distributed in the interval [— A x , A x ]. The lat- 
ter provides a simple way of implementing loss of long 
range order along the CE's. For A^, and A y we study 
the following specific cases: A x ,A y = in Sees. IIIIIIVI 
A x = A, A y = in SecEl and A^ = A v = A in 
Sers lVTIVTTI 

To study the commensurability effects, the effective 
width of the channel is varied from a value w/bo ~ 1— 10. 
We assume that the vortex density in and outside the 
channel are the same. The number of vortices in the 
channel is then given by N c h — (L/ao)(w/bo) with L 
the channel length. In a commensurate situation one has 
w = pbo and both the row spacing and (average) lon- 
gitudinal vortex spacing in the channel match with the 
vortex configuration in the CE's. When w ^ nbo these 
spacings become different, leading to generation of topo- 
logical defects. While this model differs from the experi- 
mental case where the applied field drives the incommen- 
surability, the method offers a simple way of introducing 
geometrical frustration and study various (mis)matching 
configurations. 

With a uniform transport current J applied perpen- 
dicular to the channel, the equation of motion for vortex 
i in the channel reads (in units of N/m): 

V(r) is the vortex-vortex interaction potential, j labels 
other vortices inside the channel, the damping parameter 
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7 is given by 7 = B&o/Pf with pf the flux flow resistiv- 
ity, and / = J$ is the drive along the channel. For 
films which are not too thin compared to the penetration 
depth A and magnetic fields small compared to the up- 
per critical field B C 2, the interaction V(r) is given by the 
London potential: 



V(t) = U K (\r\/X), 



(6) 



where Uq = ^q/2ttpoX 2 and <E>o is the flux quantum. 

In the simulations we integrate Eq.(|5j) numerically 
for all vortices in the channel. We use a Runge-Kutta 
method with variable time steps such that the maximum 
vortex displacement in one iteration was ao/50. Dis- 
tances were measured in units ao (? = r/ao), forces in 
units U /a and time in units ^Oq/Uq. Following Ref. 
|37| . the London potential was approximated by: 



V'(F) = In I f 



with a cut-off radius r c corresponding to r c ~ 3A/ao. 
We performed most simulations for r c = 3.33. Peri- 
odic boundary conditions in the channel direction were 
employed. For each w/bo we relaxed the system to 
the ground state for / = 0. We found that this is 
best achieved by starting from a uniformly stretched 
or compressed n or n ± 1 configuration. For an initial 
configuration with N c h vortices distributed randomly in 
the channel, relaxation resulted in (slightly) metastable 
structures, even when employing a finite temperature, 
simulated annealing method. Some peculiarities asso- 
ciated with such structures are mentioned in SecsHVl 
and IVII After the / = relaxation, the average ve- 
locity versus force (v-f) curve was recorded by stepwise 
varying the force from large f ma x - > (occasionally 
/ = — > fmax — > was used to check for hysteresis). 
At each force we measured v(f) = (xi)i t t after the tem- 
poral variations in v became < 0.5% (ignoring transients 
by discarding the data within the first 3ao). In addi- 
tion, at each force we measured several other quantities, 
e.g. the temporal evolution of ri and the time dependent 
velocity v(t) = (ai*)j. 



III. SINGLE CHAIN IN AN ORDERED 
CHANNEL 

The first relevant issue for plastic flow and commensu- 
rability effects in the channel is to understand the influ- 
ence of periodically organized vortices in the CE's (see 
Fig. The characteristic differences between com- 

mensurate and incommensurate behavior can be well un- 
derstood by focusing on a 1-D model in which only a sin- 
gle vortex chain is present in the channel. Therefore the 
CE's are assumed to be symmetric with respect to y = 
(i.e. Ax = in Fig. |3J) and only the longitudinal degrees 
of freedom of the chain are retained. At commensurabil- 
ity, w = bo, the longitudinal vortex spacing a — oq. For 



w 7^ by the average spacing a = $o/(Bw) = aobo/w does 
not match with the period ao in the edges and intersti- 
tials or vacancies develop in the channel. Their density 
c d is given by c d = |oq 1 - a _1 | = (l/a )|l - (w/b )\. 



A. Continuum sine-Gordon description 

We first consider the interaction of a vortex in the 
channel with the periodic arrays in the CE's. As shown in 
AppEl when B < 0.2B c2 and A > ao the edge potential 
arising from this interaction is: 

V ce , (x, y) = -2U e~ k ^ w+bo ^ 2 cosh(k Q y) cos k x, (8) 

where ko — 2n/aQ. For w = bo and y — 0, the associ- 
ated sinusoidal force caused by the edge has an amplitude 
which we denote by p: 



p = (47rC/ /ao)e 



[/o/(67rao). 



(9) 



Next we consider a static chain of vortices inside the 
channel. The chain is most easily described in terms 
of a continuous displacement field u(x), representing the 
deviations of vortices in the chain with respect to the 
commensurate positions, i.e. u(iao) — Ui = x% — iao- 
The edge force is then given as f p = — p sm(k u). To de- 
scribe the interaction between vortices within the chain, 
we assume that their relative displacements are small, 
d x u 1. Then one can use linear elasticity theory. Tak- 
ing into account that the interaction range A > ao, the 
elastic force at x = iao is: 



fe 



s=ja >ia 



d 2 s V(s)[u(x- 



-u(x- 



■2u(x)}. (10) 



Using the Fourier transform of V the force due to a dis- 
placement u q (x) =Ke(u q e lqx ) with wave vector q is: 



fei 



dk Uoirk 2 



2n VP + A 



_ ^ 2e lks 



[1 - cos(qs)]u(x). (11) 



Recasting this into a sum over reciprocal vectors Iko ± 
q and retaining only the / = term, one obtains the 
following dispersive elastic modulus of the chain: 



UqttX/ ap 
y/l + \ 2 q 2 



(12) 



For deformations of scale > 27rA, the elastic force is f e i = 
KqO^u with a long wavelength stiffness Ko = f^o 7r (<V a o)- 
The equation of motion for u for a uniformly driven 
chain, is obtained by adding the driving force / to the 
edge force and the intra-chain interactions resulting in: 
"fdtu = f + f p + fei- Assuming for the moment that the 
long wavelength description is valid, the evolution of u is 
given by the following sine Gordon (s-G) equation: 



jd t u = f — fj, sm(k u) + n d 2 u 



(13) 
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FIG. 4: Drawn lines: the anti-kink solution Eg. ll4H for 
X/ao = 1, 2 and 9 (most extended line). Symbols: numer- 
ically obtained displacement field for an isolated interstitial 
for the corresponding r c . Dashed line: defect shape Eg. 1181 
in the nonlocal limit. 



A useful visual representation of Eg. (jl^j) is an elastic 
string of stiffness kq with transverse coordinate u(x) in 
a tilted washboard potential (/i/fco) cos(kou) — fu. 

The s-G equation (|13[) has been thoroughly studied in 
different contexts (e.g. |38l l39l liQjl. In the static case 
(/ = 0), it has the trivial solution u — 0, corresponding 
to a commensurate chain, or kinked, incommensurate, 
solutions in which u(x) periodically jumps by ±ao, each 
jump representing a point defect in the channel. In the 
context of long Josephson junctions (LJJs, HU), a kink 
corresponds to a Josephson vortex where the phase dif- 
ference across the junction changes by 2ir. An isolated 
defect is represented by the familiar 'soliton' solution of 
the s-G model: 

Udix) = 2oo arctan[exp(±27r(a; — x c )/ld)]/ir- (14) 

Here x c denotes the center of the defect and the + (— ) 
sign denotes a vacancy or interstitial (kink or antikink). 
The length Id represent the core size of the defect: 



l d = 2iraoy/g, 



(15) 



with g the dimensionless ratio between the chain stiffness 
and maximum curvature of the pinning potential: 



g = K /2irna = 37r(A/a ), 



(16) 



as follows from Ea . . ltT2"|) . For A/a > 1, Id thus con- 
siderably exceeds the lattice spacing. The continuum ap- 
proach is validated since d x Ud < 2a$/ld <C 1. In Fig. 0]we 
have illustrated the characteristic defect shape Ea. (|14ll . 
along with numerical data from a later section. 

The long wavelength limit is only valid whe n Id c on- 
siderably exceeds A. Since Id grows only as \J X/ao the 
dispersion in the elastic interactions becomes important 
beyond a certain value of X/clq. This value is estimated 
by setting Xqd = 2ir\/ld = 1 in Ea. 1)12(1 . resulting in 
X/ao — 9, in which case Id — 54ao- 

For larger interaction range one employs the follow- 
ing approach, first derived by Gurevich 28] for mixed 
Abrikosov- Josephson vortices in grain boundaries. Ex- 
pression i|l(J|) for the elastic force can be written as an 



integral f e i = J (ds / ao)d 2 V (s)u(x + s). For defects on a 
scale < A, only the short distance behavior of V has to be 
retained: V(x) = UqK$(x/X) ~ — Uq ln(|a;/A|). Integrat- 
ing the expression for f e [ by parts and adding the edge 
force and the drive, the equation of motion becomes: 



7<9 t w = / - fi sin(fc w) + (U / a ) 



ds 



(17) 



A static solution of Eq. I|17|) for a single defect is [2£| : 

k u = 7r + arctan(±27ra/^ z ), (18) 

with Q = 6ir 2 a , which is valid when A > Q . The 
value for I T is nearly the same as the s-G core size Id for 
A/ao = 9. This means that upon approaching the nonlo- 
cal regime, the increase of core size saturates at ~ 60eto, 
while only the tails of the defect are affected according to 
Eci. (|18|) . see Fig. A more accurate calculation of the 
onset of the nonlocal field regime using Brandt's field- 
dependent vortex interaction (App^l shows that non- 
locality is only relevant for a channel in a superconductor 
with A/£ > 50. 

So far, we discussed isolated defects. For finite defect 
density, the repulsive interaction between defects of the 
same 'sign' causes a periodic superstructure in the chain. 
When Cd grows to ~ I /Id, the defects start to overlap 
significantly. For the (local) s-G model, explicit solutions 
for u have been obtained in terms of the Jacobi elliptic 
functions, for which we refer to j3^,^3|- Recently, also in 
the non-local limit where Vf > A, the 'soliton' chain has 
been described analytically , which we will not repeat 
here. 



B. Transport properties 

With a uniform drive /, the transport properties 
strongly depend on the presence and density of defects 
in the channel. At commensurability (a = do, c<j = 0), a 
threshold force f s = \i is required, above which all vor- 
tices start moving uniformly. Their velocity is identical 
to that of an overdamped particle in a sinusoidal poten- 
tial: v — yj f 2 — fi 2 |4ll [. The threshold /x coincides with 
the well known relation between shear strength and shear 
modulus of an ideal lattice by Frenkel (3(| : for a harmonic 
shear interaction, a value A = A = ao/(2nbo) = 1/7TV3 
applies in Eq.QJ. Identifying F s ao&o = fs = A* f° r 
w = bo, one finds: 



c 66 = 7rV3^/(2a ) = Ua/(8a a bo) 



(19) 



which coincides with the familiar expression for the shear 
modulus in the London limit: Cg6 = ^oB/(16tt^qX 2 ). 
In App^ we generalize the expression for the ordered 
channel potential to higher field and show that also in 
that case the potential is harmonic and that A = A° 
holds for a commensurate channel. 

At incommensurability, depinning of the chain is gov- 
erned by the threshold force to move a defect. In 
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the present continuum approach such threshold is ab- 
sent. However, taking into account the discreteness of 
the chain, in which case Ea. 1)13(1 turns into a Frenkel- 
Kontorova (FK) model 35] , a finite Pcicrls-Nabarro (PN) 
barrier exists to move a defect over one lattice spac- 
ing (see e.g. |40J). The magnitude of the PN barrier 
has been studied for a variety of cases, including FK- 
models with anharmonic and/or long-range interactions 
0,^2- For g < 1, fpN can amount to a considerable 
fraction of fi. Addition ally, in this regime anharmonic- 
ity may renormalize g |40l l42| and cause pronounced 
differences between the properties of kinks (vacancies) 
and antikinks (interstitials). In our limit g 1, where 
d x u « 1 and harmonic elastic theory applies, these 
differences are small and the pinning force vanishes as 
fpN — 327r 2 (7/iexp (— Tr 2 ^/g). Hence, defects in an or- 
dered channel give rise to an essentially vanishing plastic 
depinning current J s [43T |. 

Considering the dynamics, for small drive / < fx the 
motion of defects, each carrying a flux quantum, provides 
the flux transport through the channel. When defects are 
well separated, for Cd < l^ 1 , the mobility of the chain is 
drastically reduced compared to free flux flow and the 
average velocity v is proportional to the defect density: 
v = CdVdao. Here Vd is the velocity of an isolated defect 
at small drive. It can be calculated from the general 
requirement that the input power must equal the average 
dissipation rate: 



/ = jv + (rfv/ao) / (dshyds 



(23) 



fv = i{{d t u) 2 ) LA = ( 7 /i d ) / {d t u) 2 dx 



(20) 



The last step arises from the space and time periodicity 
of u. Using dtu — Vdd x u and the kinkshape Ea.(jl4p. one 
obtains the 'flux flow resistivity' at small defect density: 

dv/df = c d ao(ir 2 ^g-/2 1 ) 7 (21) 

where ir 2 ^/g/2-f = Md is the kink mobility in the s-G 
model |3j| • For larger defect density, where defects start 
to overlap this relation changes. The linear response for 
/ < fi may then be obtained from the solutions for u 
based on elliptic integrals |28l 1381 |40j . 

For larger drive / > fi, the 'tilt' induced reduction of 
the (washboard) edge potential becomes important. This 
leads to an expansion of the cores of the sliding defects 
and causes a nonlinear upturn in the v-f curves. Exact 
solutions of EqE3 describing this behavior do not exist. 
Therefore we use a perturbative method similar to that 
in [2^, 0| which is able to describe the full v-f curve 
over a wide range of defect densities. It is convenient to 
define the displacements h(x, t) = u(x, t) — s(x, i), where 
s(x, t) = (q/ko)x + vt, with (q/ko) = c^ao, is the continu- 
ous field describing the displacements of an undeformed 
incommensurate chain (i.e. straight misoriented string 
in the washboard potential) moving with velocity v. In 
terms of h, the equation of motion (|13|) and Eq. 120|) can 
be written as: 

7«(1 + d s h) = f + fisin(k h + qx + k vt) + n(q/k a ) 2 d 2 h 

(22) 



The last term in Eq. I)23|l describes additional dissipation 
due to internal degrees of freedom in the chain. Under 
influence of the potential, h acquires modulations with 
period 1/cd in x, i.e. period ao in s. These modulations 
are then expressed as a Fourier series of modes with wave- 
length l/(mc d ) (m integer > 1) and amplitude h m : 



h(x, t) = h, 



, exp[jm/cos] + c.c. 



(24) 



The overlap of defects and the core expansion for / > /i 
appears in the q and v dependence of h. Both effects 
cause a reduction of the relative displacements h. An 
approximate solution for h(v) is obtained by substituting 
Ea. (|24|) into Eq. l|13fl . yielding the coefficients h m (details 
of the solution are deferred to App lB"Jl . The v-f relation 
Ea. (|23|) attains the form: 



/ = jv 



2M? 



(25) 



The additional 'friction' force is represented in terms of 
the pinning frequency lu p = the washboard fre- 

quency ujq = kov and uj r — K 2 ^{cd)/^ which is the effec- 
tive relaxation frequency for nonlinear deformations as- 
sociated with a defect density Cd = qj (27r), with K 2 ^(cd) 
given in App[5] At small v, the elastic relaxation time 
\/uj r for the chain to relax is much smaller than the 
timescale l/wo between passage of maxima in the edge 
potential. This corresponds to the linear sliding response 
of the static structure of (overlapping) defects. For large 
v , l/uj r 3> l/^o meaning that the incommensurate chain 
is not given enough time to deform. This leads to ex- 
panded defects described by a sinusoidal variation of h 
with reduced amplitude (see App[5J|. The v-f curve then 
approaches free flux flow according to / — 71; ~ w _1 as 
for a single particle. 

Recently, exact solutions describing the nonlinear 
dynamics and core expansion of mixed Abrikosov- 
Josephson vortices based on the nonlocal Ea. l|17[l have 
been derived in |28j . The resulting transport curves are 
very similar to those obtained from Ea. (|25|l . see Fig. [S] 
We also note the similarity with the IV curves obtained 
from a model for kinked Josepson strings |45| in high 
T c superconductors with the field under an angle with 
respect to the insulating layers. 



C. Numerical results 

The simulations of symmetric channels (Ax = 0) for 
w ~ bo fully support the above findings. The interaction 
with the CE's for w — bo provides a maximum restoring 
force with a value 0.054, independent of the interaction 
cut-off r c used in the numerics. This value is in agreement 
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FIG. 5: v-f characteristics for ordered vortex channels with 
w w foo and Ax = 0. Symbols are simulation results, drawn 
lines are obtained with Eg. 1251 . The inset shows an expanded 
view of the small velocity regime. 




FIG. 6: Displacements u(x) along the channel for w = 
0.97bo:(°) Numerical result for / = also representing a 
snapshot of the moving chain at low drive, / = 0.01. The 
thick drawn line shows the result for u as calculated from 
the Fourier modes given in App|B| (•) Displacement field for 
v = 0.09 (/ = 0.1). The data mask a drawn line which is 
obtained from Eq. lB8t in App|B| The straight drawn line 
shows the displacement field s(x) in absence of the periodic 
potential. 



with the dimensionless values for \x and cee in Eq. © and 
Cd: aofi/U = 1/(6tt) and c 66 a 2 /U = 1/(4^3)). 

The data points in Fig. 0]show the displacement field 
of a single defect (obtained by adding one vortex to a 
commensurate chain) for three values of the cut-off r c 
(i.e. various A/a ). We conclude that up to r c = 30 
(A/ao = 9) the s-G kink shape Ea. lpH|) forms a good 
description of a defect in the chain. 

The data points in Fig. [3] show numerical results 
for the transport of a single chain in channels of vari- 
ous widths and r c = 3.33. The features discussed pre- 
viously, i.e. the vanishing PN barrier and nonlinear 
transport, clearly appear in the data for incommensu- 
rate chains. We also plotted the results according to 
Ea. (|25|l . with Rifled) evaluated using the results in 
AppIbI for A/ao = 1 and taking into account that \x 
slightly depends on w. The analytical treatment gives 
a very reasonable description of the data. Finally, we 
show in Fig. the numerical results and analytical re- 
sults of App[5]for the quasi-static and dynamic shape of 
the chain for w/b — 0.97 (q = 0.03/a ). The numerical 
results closely mimic the analytic results, both for the 
kinked shape at small v and the core expansion with the 
associated reduction of h for large v. 

To conclude this section we mention that, at incom- 
mensurability, due to the vanishing barrier for defect 
motion, the average velocity (ii(t))i has a vanishing ac- 
component. Only at commensurability the washboard 
modulation is retained, the velocity at large drive being 
(v(t)) = v + (fj,/-y)sm(u) t). 



IV. ORDERED CE'S AND MULTIPLE CHAINS 

We now turn to the results for channels containing 
multiple vortex rows and ordered CE's. The simulations 
are performed with the full 2D degrees of freedom and 




f(U /a ) 

FIG. 7: Transport curves for commensurate channels with 
w/bo = 2, 3, 4, 5, 7, 9 from right to left. The thick drawn line 
represent free flux flow. Inset: friction force / — 'yv vs. v 
for w/bo = 2, 3, 4, 5, 7, 9 from top to bottom. The dotted line 
represents Ea, l28H for n — 9. 

r c = 3.33. We implemented an edge shift Ax(w) with a 
saw tooth shape (0 < Ax < ao/2). This assures that, as 
we vary in, a perfect hexagonal structure is retained for 
w = pbo with p an integer. However, for w ^ pbg, the 
qualitative behavior did not depend on Ax. 

Figure \7\ shows v-f curves of commensurate channels, 
w/bo = n with integer n > 2. In these cases the ar- 
rays are perfectly crystalline and have a shear strength 
f s = fibo/w, inversely proportional to the channel width 
and in accord with Eq.QJ with A = A = This is 

consistent with the fact that only the first mobile chains 
within a distance ~ bo from both CE's experience the pe- 
riodic edge potential (see Eq. (|A3|l ) while the other chains 
provide an additional pulling force via the elastic interac- 
tion. This interaction brings an additional feature to the 
dynamics, namely shear waves. The shear displacements 



of rows n in the bulk of the channel can be described in 
continuum form, u n (t) — > u(y,t), by the following equa- 
tion of motion: 



w= 3.92b n 



jd t u(y,t) = f + c e6 a badyu(y,t). 



(26) 



At large v the CE interaction can be represented by os- 
cillating boundary conditions. As shown in AdpIUI this 
causes an oscillatory velocity component dh/dt with y- 
dependent amplitude and phase describing periodic lag- 
ging or advancing of chains with respect to each other: 

d t h(y, t) f(y) sm(uj t) - g(y) cos(u t). (27) 

Here loq is the washboard frequency kov, f(y) = 
cos(y/l±, v ) cosh(y/l± lV ), g( y) = sm (y/l± iV ) sinh(y/Zj_ )t) ). 
The length scale Zj_,u = yj (fi/jv)^ explicitly depends 
on v and represents the distance over which the ampli- 
tude and phase difference decay away from the CE's. Al- 
though in principle Ea. i|27|) is only valid for jv/fi > 0.25, 
it provides useful qualitative insight in the dynamics at 
all velocities: at small velocity l± iV is large, meaning that 
for all rows the velocity modulation and phase become 
similar. Hence, for v — ► the array may be described as 
a single vortex chain, which is the underlying origin of 
the fact that close to threshold the curves approach the 
ID commensurate behavior v = \J f 2 — ff with reduced 
threshold f s = /i/n. At large velocity, Zj_ jt , eventually 
becomes less than the row spacing. In that limit only 
the two chains closest to the CE experience a significant 
modulation. In AppOwe quantitatively analyze the fric- 
tion force in this regime with the result: 



f -yu = 



(28) 



In the inset of Fig. [71 this behavior is displayed for n = 9 
by the dotted line. In the high velocity regime the result 
agrees well with the numerical data, at lower velocities 
Ea. (|28|l underestimates the true friction. 

Next we discuss the behavior of incommensurate chan- 
nels. The static vortex configuration for a channel of 
width w/bo = 3.92 is shown in the upper part of Fig. [S] 
A Delauney triangulation shows that the array consists 
of 4 rows with two pairs of 5, 7- fold coordinated vortices 
at the CE constituting two misfit dislocations of opposite 
Burgers vector b and glide planes along x. Due to their 
mutual attraction, dislocations at the upper and lower 
CE are situated along a line with an angle of ~ 60° with 
x. The two edge dislocations thus form a 'transverse' 
stacking fault (TSF). In the lower part of Fig. the 
structure for a channel with w/bo = 3.52 is shown. Here 
the density of stacks, given by ctsf — |(V a o) — (V a )|) 
is enhanced. The dislocations at one side of the CE re- 
pel each other and are equally spaced, like the periodic 
superstructure for a single chain in Fig. The slight 
misalignement between the 'upper' and 'lower' disloca- 
tions of a pair is due to the relative shift between the 
CE's: the exact orientation of the pairs is determined 




FIG. 8: Delauney triangulation of the static structure for 
two incommensurate channels: w/bo = 3.92 and w/bo = 3.52. 
Open circles and filled squares denote 7 and 5-fold coordinated 
points, respectively. The construction for the Burgers vector 
is shown for w/bo = 3.92; the drawn lines for w/bo = 3.52 
mark the transverse stacking faults. 
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FIG. 9: v-f curves for incommensurate channels. Inset: 
regime of small drive. 



by the choice of Ax. For channel widths in the regime 
n < w/bo < n + 1/2 with integer n, we find very simi- 
lar structures but instead of TSF's consisting of vacan- 
cies, we now have TSF's consisting of interstitial vortices, 
again arranged in a periodic superstructure. 

Figure shows the transport curves associated with 
these structures. As for the single chain, the presence 
of misfit defects causes an essentially vanishing threshold 
force. For small drive, / < [ibo/w, a low mobility regime 
occurs associated with glide of the edge dislocation pairs 
along the CE. This allows for elastic motion of a complete 
TSF, i.e. the vortices in the 'bulk' of the channel remain 
6-fold coordinated. 

It is interesting to study how the mobility due to 
the TSF's changes on increasing the number of rows. 
In Fig. EH we plot the mobility per stack, Mtsf — 
(dv/df) f^o/cTSF versus channel width. Mtsf around 
each peak decreases with increasing ctsf- This is caused 
by overlap of the strain fields of the defects, in analogy 
with the behavior for a single chain. The overall increase 
of the peak value of Mtsf is related to a change in the 
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FIG. 10: The mobility per stack Mtsf = {dv / df) / c stac k ver- 
sus w/bo. The drawn line shows the predicted form Ea. l|29|l 
for A/ao = 1. 




FIG. 11: Threshold force versus w/bo for ordered channels. 
The dashed line represents Frenkels prediction for an ideal 
lattice in the continuum limit. 



size of an isolated TSF. An extension of the analysis 
in Sec lIIII allows to describe this change quantitatively. 
For small n the longitudinal deformations do not vary 
strongly over the channel width. This can be understood 
by considering shear and compression deformations re- 
lated by the equation nd^u x + ceeaobodyU x — 0. It fol- 
lows that a longitudinal deformation on a scale h\ along 

the channel varies over a scale l± = hi y/ ceeaobo/K per- 
pendicular to the channel. In case l± > w, the transverse 
variation of u x (y) is small and can be neglected so that 
kq in Eq. (|13|1 can be replaced by an effective stiffness ukq 
due to n rows. Similarly, the driving force is replaced by 
/ — > nf. This results in the same equation (|13|) with a 
rescaled edge force /i — > fi/n. Accordingly, the longitu- 
dinal size of a defect (TSF) is given by Itsf = Z^ao^/rig 
and the mobility of an isolated TSF by (compare Md 
below Eq.lEU)): 

M TSF tt 2 ^/2 1 . (29) 

As shown by the drawn line in Fig. 1101 this form gives a 
reasonable description of the data up to n — 3. Working 
out the condition l± > w given above for the validity 
of Eq.|J55J, one obtains w < (ld/2)\/c e6 aobo/ k — 3&0) in 
agreement with the data. At larger n, 'bulk-mediated' 
elasticity [3|| leads to decay of the longitudinal deforma- 
tions towards the channel center. We also note that, due 
to the increase of Itsf with n, the density ctsf for which 
defects are non-overlapping, decreases on increasing n. 

In the v-f curves of Fig. Elwe observe at larger velocity 
features very similar to the transport of the ID chain: 
f° r / > aV 71 the effective barrier is reduced, leading to 
core expansion of the TSF's. Accordingly, the curves 
approach free flux flow behavior. As in the commensurate 
case, this approach is initially more slowly then / — 'yv ~ 
1/v due to additional oscillating shear deformations in 
the channel for / > fi/n. 

Figure 1111 summarizes the behavior of the shear force 
f s , taken at a velocity criterion v ~ O.OI/1/7, versus the 
matching parameter. At integer w/bo — n, the threshold 
is f s = n/n, but we note that it can be reduced due to a 
finite edge shift Ax. At mismatch f s is essentially vanish- 
ing. Near 'half filling', w/bo ~ n± 1/2, where the arrays 



switch from n to nil chains, a small enhancement of f s is 
observed. In this regime, the static (/ = 0) structure was 
obtained by annealing from a random initial configura- 
tion, attempting to determine the exact switching point. 
This results in metastable structures with coexisting n 
and nil row regions (or longitudinal stacking faults, 
LSF's) bordered by dislocations with misoriented Burg- 
ers vector, see 35]. The increase in f s is caused by the 
finite barrier for climb-like motion of these dislocations, 
via which an LSF can move as 'giant' defect through the 
channel. For sufficiently large drive, (part of) the LSF's 
are annealed which may result in hysteresis for up/down 
cycled v-f curves. We will discuss these 'mixed' n/n i 1 
structures in more detail in Secs lVIl and IVIII in the con- 
text of disordered CE's. We also note that the integer 
chain structures with TSF's away from half filling differ 
from the results in |35j. The structures there, obtained 
from a random initial configuration, contained point de- 
fects unequally distributed among rows, yielding 'gliding' 
dislocations within the channel. Such structures are also 
slightly metastable but the conclusion of vanishing /■ for 
incommensurate, integer chain structures, drawn in [3J| . 
remains unaltered. 



V. SINGLE CHAIN IN A DISORDERED 
CHANNEL 



We will now consider the influence of disorder in the 
CE arrays on transport in the channels, focussing in 
this section on the characteristics of a single chain for 
w/bo ~ 1 with only longitudinal degrees of freedom. The 
CE disorder is implemented with longitudinal random 
shifts as described in Sec[H] We note that both CE's 
remain 'in phase'; the effect of quenched phase slips or 
dislocations between the CE's will be treated in the dis- 
cussion in Sec lVIIll 
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A. Disordered sine-Gordon equation 



coupling to the strain: 



First we consider the form of the channel potential in 
presence of weak disorder. To that end we generalize 
Ea. (|Al|l in ApplAland express the CE potential at ro = 
[x, y — 0) in terms of the vortex density p e in the CE's: 



Ke(ro) = (2tt)- 2 / aW(k)p e (k)e lkr °, 



(30) 



with p e (k) the Fourier transform of p e . For weak dis- 
order (V ■ d <C 1), this density can be expressed in 
terms of the displacement field d in the CE as fol- 
lows H3: p e (r e ,d) ~ (B/® )(1 - V • d + Sp e ) where 
Sp e = ^cos[Kj(r e — d(r e ))] represents the microscopic 
modulation due to the lattice (K^ spans the reciprocal 
lattice) while V • d reflects density modulations. As de- 
scribed in AppIdI this decomposition of p e leads to two 
contributions to the potential: 

V ce = Vi(x)+V p (x) 

= -(B/$o) J dr e V(r - r e )V • d(r e ) 

-\p + 5p(x)] cos[fc (a; - d)]/k , (31) 

where in the second term Sp(x)/p = iT\/Sd x d. The 
term Vi represents long range potential fluctuations and 
is smooth on the scale ~ a . Its correlator T;(s) = 
(Vi(x)Vi(x + s)} is derived in AppJDJ Assuming that d x d 
has short range correlations (on the scale ~ cto/2) and a 
variance ((d x d) 2 ) = A 2 /3 as in the simulations, Ti can 
be written as: 



r,(a) ~ C a A 2 U 2 (X/a ) 



ll+«p-(x) 2 



(32) 



The exponent a and the prefactor C a depend on the 
disorder correlations between rows in the edge: a — 2 
when the strain d x d{x) is identical for all rows and a 
1 when the strain is uncorrelated between rows. The 
term V p in Ea. l|31|) is the quasi-periodic potential arising 
from Sp e of the vortex rows nearest to the CE's. The 
amplitude fluctuations dp/ko are characterized by (see 
App|Dj: 



r«(«) 



(Sp(x)Sp(x + s)) x 



k 2 



(pAa /2) 2 e~ i ^ ) (33) 



To obtain the energy of the vortex chain and the equa- 
tion of motion, the vortex density inside the channel, p c , 



is decomposed similar to p e : p c (x,u) 



L [l -d x u 



Sp c (x, u)], where u is the displacement field of the chain. 
As shown in App[Dl in the limit A > ao the result- 
ing interaction with the CE's can be written as H = 
H SG + H a + H s where H SG = a^ 1 J dx[( Ko /2)(d x u) 2 - 
(p / k Q ) cos(k u)] represents the sine-Gordon functional 
for an ordered channel, and H a , H s are the disorder 
contributions due to amplitude fluctuations and random 



dx Sp(x 

ao fco 
dx 



■ cos(fco") 
V s (x)d x u. 



(34) 



The term V s (x) = Vi(x) — Kod x d(x) contains contri- 
butions from local and nonlocal strains. The latter 
dominates for A > ao (see App ID|) . Hence T s (s) = 

(v s (x)v s (x + S )) ~r,(«). 

The model described by H — Hsg + H a + H s is also 
used to describe LJJs or commensurate CDW's with 
weak disorder, however with different disorder correla- 
tions. In the former case, the term H a in Ea. (|34|) de- 
scribes local variations in the junctions critical current 
|48L fljil ] . For CDW's, a disorder contribution of the form 
H a arises from so called backward scattering impurities, 
while the term H s originates from 'forward' scattering 
impurities [5(J • We also note that our model differs from 
the usual Fukuyama-Lee-Rice model for CDW's [HI]], in 
which commensurability is ignored either due to strong 
direct random coupling to u (8p(x) 3> p) or due to large 
mismatch. 

In principle, the equation of motion for the chain is 
given by "fdtU = —SH/Su. However, it has been shown in 
previous studies 0, EH 01 that in the moving state 
a convective term —jvd x u should be included. While 
irrelevant for the depinning process, such term can be 
important for the dynamics and for completeness we in- 
clude it [55| . The resulting equation of motion is: 



"fd t u = f + n a d 2 u — [p + Sp(x)] sm(kou) 
-d x V s - jvd x u. 



(35) 



In writing Eq. (|35|l we have assumed, for simplicity, that 
the elastic deformations in presence of disorder can be 
described by the long wavelength stiffness ko- Ignoring 
the last term, Eq. l|35|l describes the transverse displace- 
ments u(x) of an elastic string in a tilted 'washboard' 
potential with random amplitude p(x)/ko and random 
phase 4>(x) — — f* dx'V s {x')/ n . The latter represents 
a u independent random deformation of the chain. 



B. Numerical results 

The influence of disorder on the threshold force and 
the dynamics of the chain are directly visible in numerical 
simulations. The simulations were performed using r c — 
3.33ao and channels of length L > lOOOao. 

Figure 1121 shows several v-f curves for channels with 
0.93 < w/b < 1.1 at a disorder strength A = 0.025. 
We first focus on the result for the commensurate case 
w = b{). The disorder leads to a significant reduction of 
the threshold f s with respect to the pure value p = 0.054. 
The reduction is enhanced on increasing A, as shown 
for A = 0.05 in the inset. The origin of the reduction 
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FIG. 12: v-f curves of a commensurate and various incom- 
mensurate chains for weak disorder A = 0.025. Inset: v-f 
curve of a commensurate channel for A = 0.05 



is that disorder lowers the energy barrier for formation 
of vacancy/interstitial (kink/antikink) pairs in the chain. 
FigurelT3*k shows the time evolution of the displacements 
Ui = Xi — iao upon a sudden increase of / to a value 
/ = 0.049 > f s at fx. For t < h u is 'flat' and the 2D 
crystal formed by the chain and the CE's is topologically 
ordered. At t = t\, the motion starts at an unstable site 
(at x/ao ~ 500) by nucleation of a vacancy/interstitial 
pair visible as steps of ±ao in u. We henceforth denote 
the force at which this local nucleation occurs by /„. The 
defects are driven apart by the applied force and when 
their spacing becomes ~ l d a new pair nucleates at the 
same site. This process occurs periodically with rate R ni 
leading to the formation of a domain with defect density 
c d = R n /(v d ) and a net velocity v = c d (v d )a = R n a 
with (vd) the average defect velocity. In the present case 
of weak disorder (v d ) is essentially the same as for A = 0. 
For a further increase of the force to f — 0.053 an in- 
crease of the nucleation rate is observed. In Ref.[33j we 
showed that in larger systems coarsening occurs in the 
initial stage of depinning due to a distribution of un- 
stable sites. However, after sufficiently long times the 
stationary state consists of one domain around the site 
with the largest nucleation rate (smallest local threshold 
/™ m ) with vacancies travelling to the left and intersti- 
tials to the right. It is interesting to compare this to a 
study of CDW's with competing disorder and commen- 
surability pinning [52j . U sing a coarse grained version of 
Ea. (|35[l it was found in |52j that in the pinning domi- 
nated, low velocity regime, the so-called interface width 
W(L) = \J ((u(x) — (u)) 2 ) x grows linearly with the sys- 
tem size L. The mechanism of defect nucleation which 
we observe naturally explains this phenomenon. In ad- 
dition, we found that at depinning the average velocity 
v = R™ ax a can be described by R™ ax oc (/ - f^ in f 
with a depinning exponent /3 = 0.46 ± 0.04, similar as 
previously reported for ID periodic media |56j . 

The defective flow profile does not persist up to a arbi- 
trary large forces. In the commensurate v-f curve in Fig. 
E3and its inset, a small kink is observed for / < //. Asso- 




FIG. 13: (a) Evolution of longitudinal displacements Ui(t) for 
the commensurate chain in Fig. I12f A = 0.025), plotted for 
clarity in a transverse way vs. x. At t = ti the force is 
increased above threshold, (b) Stationary evolution of u;(i) 
for large drive (/ = 0.08) showing the motion over a distance 
~ 1.8ao (the t and u axis have arbitrary offset, and a few 
frames around t = 10 were omitted for clarity). 



ciated with this kink we find a transition to a much more 
ordered state. We have illustrated the temporal evolu- 
tion of vortex displacements in this state in Fig. I13f b) 
for / = 0.08. The 'staircase' structure has vanished and 
the relative vortex displacements are greatly reduced. In 
fact, in the above mentioned study Ref. [52j a very simi- 
lar transition in the CDW dynamics was found, and was 
shown to be of first order. We leave the precise depen- 
dence of this transition on disorder and vortex interac- 
tions in our channels for future studies. 

We now turn to the incommensurate case. The v-f 
curves with w ^ bg in Fig. ^| all exhibit a finite thresh- 
old instead of the vanishing threshold for the incommen- 
surate channels without disorder (Fig. [SJ. With disorder 
the defects that are present in the channel, couple to the 
disorder in the CE, which causes a pinning barrier f d . 
This barrier has a distribution along the channel {f d } and 
maximum value f r d ax . We now focus on the curves with 
small defect density c d < l/l d for which the defects are 
individually pinned. In this regime the threshold force 
f s satisfies f s < f2 iax - 

The precise threshold behavior depends on the distri- 
bution of barriers {f d }, similar as for LJJ and CDW sys- 
tems 01- A s an illustration, we show in Fig. I14f a') the 
evolution of displacements for a channel with w/bo = 0.99 
for a force just above threshold f(t > 0) > f s . The static 
configuration at t = (/ < / s ) shows that the disorder 
breaks the periodicity of the 'soliton' chain. For t > 0, de- 
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FIG. 14: (a) Evolution of displacements for w/bo = 0.99 and 
/ = 0.013 > f a for t > 0. Defects at x ~ 50ao, x ~ 150ao and 
x ~ 325ao are initially pinned while the others are mobile. 
A defect collision-release process occurs at x ~ 150ao and 
t ~ 1500. (b)Evolution of displacements when / is suddenly 
increased to / = 0.048. Nucleation is observed for x « 500ao. 



pinning starts with the defect at x ~ 270ao and proceeds 
via a 'collision-release' process between the moving de- 
fect and its pinned neighbor. Thus, for / > / s strong lo- 
cal variations in the defect mobility exist and the overall 
chain velocity depends strongly on the distribution {fd} 
(of which we show an example below). However, as seen 
in the v-f curves in Fig. ^] for / > 2f' d nax these effects 
vanish and the mobility approaches dv/df ~ Cda^Md with 
Md the defect mobility without disorder. Another feature 
in the v-f curves for small defect densities is the velocity 
upturn at a force / ~ /™ m < M- It is caused by nucle- 
ation of new defect pairs in the incommensurate chain. 
The start of such a process is illustrated in Fig. I14f b): 
at x ~ 500ao the chain is unstable against pair nucleation 
and the nucleated interstitials/ vacancies are formed 'on 
top of the moving incommensurate structure. This pro- 
cess only occurs at small defect densities when the time 
between passage of existing defects ~ \/{vdCd) exceeds 
the nucleation time R^ 1 - For / > /i, the structure of 
both defects disappears again. The resulting dynamic 
state resembles that of the large velocity profile shown in 
Fig. El but with additional 'roughness' due to the weak 
CE disorder. 

The v-f curves in Fig. ^]at w/bo = 0.93 and w/bo 
.1.1. for which the defect density in the chain is larger, 
exhibit a smaller threshold force. In this regime the in- 




FIG. 15: Threshold f 3 , obtained from a velocity criterion v w 
0.025^i/7 and L = 1000ao, versus w/bo for several disorder 
strengths. Data were averaged over 5 disorder realizations. 
Inset: v-f curves in the strong disorder regime A = 0.2. 



teraction between defects starts to become important and 
f s is determined by collective pinning of the defects. This 
situation was studied analytically for the case of Joseph- 
son vortices in a disordered LJJ in [43. We will not 
consider this situation explicitly but we note here that, 
as the disorder and the typical pinning force on the de- 
fects increases, the onset of the collective pinning regime 
shifts to larger defect density, where defect interactions 
are stronger |49| . 

In Fig. E| we show the dependence of f s on channel 
width, both for the weak disorder regime treated above 
and for larger disorder. The data at A = 0.025 exhibit a 
sharp peak at w — bo, reflecting the gap between mini- 
mum nucleation threshold and maximum defect pinning 
force. Larger disorder however rapidly smears the peak, 
being eventually completely suppressed for A > 0.15. 
The origin of this behavior is a spontaneous nucleation 
of defects in the static chain at larger disorder. This is 
conveniently illustrated via the changes in the distribu- 
tion of the individual defect pinning force {fd} and that 
of the nucleation threshold {f n } for increasing disorder 
strength, shown in Fig. The data were obtained by 
simulating hundreds of short channels (L — 100ao) both 
with one vacancy, yielding {fd}, and without 'geometric' 
defects, yielding {/„}. While for A = 0.025 the distri- 
butions are separated (formally, in infinite systems, such 
separation only exist for bounded disorder, see the next 
section), for larger disorder they start to overlap and they 
become nearly identical for A = 0.075. This implies that 
nucleated defect pairs at w — bo can remain pinned, while 
at incommensurability defects may be nucleated before 
the 'geometrical' defects are released. In other words, 
regardless of the matching condition the static configu- 
ration contains both kinks and antikinks. 

While for bounded disorder static defects first appear 
at a disorder strength defined as A c , the complete col- 
lapse of the peak in f s is associated with the presence of 
a finite density of disorder induced defects, of the order 
of the inverse kinkwidth ~ 17 . We define the disorder 
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FIG. 16: Probability density of critical forces in channels of 
length L = lOOao ~ 5ld and w — bo for a commensurate chain 
(•) and for a chain with one vacancy (o) for A = 0.025 (a), 
A = 0.05 (b) and A = 0.075 (c). 
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FIG. 17: Lower panel: relative displacements m — df for A = 
0.15 and / < f s — 0.0165 at commensurability w — bo. Upper 
panel: evolution of displacements at depinning, / = 0.017. 
The time increment between consecutive snapshots is At = 10 
and for clarity each snapshot has been shifted up by ao- 



strength at which this occurs as A*, here A* ~ 0.15. 
An example of the displacement fields for this disorder 
strength is shown in Fig. El f° r w — bo. The lower 
panel shows the displacements for / < / s , relative to the 
displacements in the CE. Clearly, the static configura- 
tion has numerous defects. In general, the approach to 
the critical pinned state occurs by avalanches in which 
local nucleation and repinning, i.e. nonpersisting nucle- 
ation events, drive the rearrangements. The upper panel 
displays the evolution of displacements above threshold, 
revealing a growth of 'mountains' due to persistent nu- 
cleation, superimposed on a disordered background. 

The effect of large disorder on the shape of the v-f 
curves is shown in the inset to Fig. 1151 All curves now 
exhibit essentially linear behavior |57j . except in a small 



regime just above f s . We note that also in this disordered 
regime a gradual transition to a smoother displacement 
field occurs at larger forces, similar to the dynamic tran- 
sition found for CDW's. Going back to Fig. EDwe should 
also mention the overall asymmetry of f s with respect to 
w/bo = 1 and the slight decrease of f s on increasing w/bo 
at large disorder. These effects are unrelated to the com- 
petition between commensurability pinning and disorder 
discussed so far, but simply reflect the overal decrease of 
the edge potential for larger width. 



C. Analysis of pinning forces and crossover to 
strong disorder 

Using the results of Sec IV Al we now analyze in more 
detail the dependence of the pinning force on disorder 
and the vortex interaction range. We focus on the average 
pinning strength of isolated defects, which we derive here 
in a semi-quantitative fashion (the formal calculation is 
deferred to App[DJ). Our analysis applies to the case 
of weak disorder, i.e. we assume that the defect shape 
is unaffected by disorder H^)- Extrapolation to larger 
disorder provides a useful estimate for the crossover value 
A* at which the commensurability peak vanishes. We 
conclude the section with a summary of previous results 
[33| for the threshold forces /™ n and f™ ax in the special 
case of bounded disorder. 

The disorder correction Ea. (|34|l to the energy of the 
vortex chain consists of a term H a , due to amplitude 
fluctuations in the periodic potential, and a term H s , due 
to random coupling to the strain. We fir st evaluate the 
typical pinning energy of a defect \J (E 2 ) due to the am- 
plitude fluctuations. The local fluctuations are assumed 
to be uncorrelated on a length scale ao, and have a vari- 
ance ((5 /j, / ko) 2 } . Hence, for a defect in the chain, which 
extends over a range l d , the resulting random potential 
has a variance (E%) ~ ((5/x/^o) 2 )(W a o) — [i 2 lS?l d ao/ ^- 
The typical pinning force on a defect is then given by 
\J {E 2 ) /l d which reduces to: 



y(/f>^0.2 M A. o - 1 /4 



(36) 



The typical pinning energy \J {E 2 ) of a defect due 
to the term H s in Eq.(j2U is estimated in a similar 
way. The mean square energy due to coupling of a 
single fluctuation in V s to the strain of a defect is 
~ T s (0)(r d /a ) 2 (2a /l d ) 2 where r d is the range of T s , 
given below Ea.l|M)l. and a /l d represents the strain. On 
the scale of a defect, there are l d /r d such fluctuations. 
Thus the associated random potential for a defect has 
a variance {E 2 ) ~ T s (0)(r d /a o ) 2 (a o /l d ) 2 (l d /r d ). Taking 
r s (0) ~ T;(0) and using Ea. (|32J) . in which case r d = A, 
yields (E 2 ) ~ 2C a (U AX) 2 (X/a ) a /{l d a a ) (the factor 2 
comes from the refined calculation in App[DJ The typ- 
ical pinning force y/ (E 2 ) / l d due to random coupling to 
the strain is then given by: 



(37) 
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The ra tio b et ween t he two characteristic defect ener- 
gies is 

V / MI/v / My - 8 v / C i I(A/ao) H * i , which shows 
that, particularly for increasing A/ao, the dominant pin- 
ning is due to random coupling to the strain. Hence- 
forth we use only this contribution. We next estimate 
the disorder strength where the commensurability peak 
vanishes. As mentioned before, this collapse occurs when 
the density of disorder induced defects becomes ~ Zj 1 , 
in other words, when the typical energy gain of a de- 
fect due to disorder becomes of the same magnitude as 
its bare elastic energy J(dx/a )(K /2)(d x u) 2 ~ fia 0y /g. 

This leads to A, oc 1/2 cT (2q+1)/4 . For the particu- 
lar case of random strains that are identical for all rows 
(a = 2), A* is given by: 

A, ~ 3g- 5/4 (38) 

This can be compared to the numerical data in Fig. IT51 
Even though those results were obtained for A/ao = 1 
(g = 3tt), formally outside the regime of validity of our 
analysis, the predicted value A* ~ 0.18 is in reasonable 
agreement with the data. 

In the particular case of bounded random strains in 
the CE's, the distribution of the nucleation force {f n } at 
commensurability is bounded from below by /™ m and 
that of the defect pinning force {fd} is bounded from 
above by /™ a:E (at weak disorder). For completeness we 
give here the previously derived results 33] for these ex- 
tremal values: both occur due to disorder fluctuations on 
the same length scale as that on which the displacement 
field u(x) varies. For a defect this naturally corresponds 
to The associated maximum defect pinning force is 
f™ ax I V oc Agr 3 / 2 (for uniform strains |23). For nucle- 
ation, at / < (j,, the appropriate length scale is l san , the 
extent of a so called small amplitude nucleus [3^ ■ Due to 
the nonlinearity of the pinning force, l san itself depends 
on the force, i.e. I sm (f]> Id and it diverges for / — > //. 
As shown in detail in [33j, this leads to a minimum nu- 
cleation threshold given by 1 — (/™ to /m) oc [c/ 3 / 2 A] 4 / 3 . 
From the condition f d nax = /™ m one then obtains the 
disorder strength A c at which pinned defects can first 
appear spontaneously in the system: A c ~ 5~ 3 ^ 2 < A*. 



VI. WIDE CHANNELS WITH WEAK 
DISORDER 

We now consider how channels of larger width, in 
which vortices have the 2D degrees of freedom, behave 
in the presence of weak edge disorder. Close to com- 
mensurability, the effects we find are similar to that for 
a single chain. However, around 'half filling' the impor- 
tance of the transverse degrees of freedom of the channel 
vortices become apparent. 




FIG. 18: v-f curves for weak disorder (A = 0.05) and several 
channel widths. The arrow indicates the yield strength f s = 
j«/3 for w/bo = 3 and no disorder. All data are computed for 
channel lengths L > 400ao. 



A. Behavior near commensurability 

For the commensurate case, w/bo — n, weak CE dis- 
order causes a reduction of the threshold with respect to 
the ideal value /° = /u/n, see the data in Fig. ED The re- 
duction originates from defect formation at threshold, as 
illustrated for w = 3&o in Fig. El In (a), three snapshots 
of the displacements of individual rows inside the channel 
are displayed. The first snapshot is for / < / s ~ 0.7/° 
and yields the 'flat' profile. The subsequent snapshots 
for f > f s reveal simultaneous nucleation and motion of 
a pair of 'oppositely charged' TSF's, each terminated at 
the CE's by a pair of edge dislocations (see Fig. Hffl b)). 
The macroscopic, stationary motion of the array is gov- 
erned by periodic repetition of this process at the least 
stable nucleation site. In (c) we show the vortex trajec- 
tories during nucleation of the TSF at the left in Figs. 
I19f a).fb'). Very similar images were obtained in deco- 
ration experiments at the initial stage of VL depinning 
in NbSe2 implying that even for weakly disordered 
VL's, defects may nucleate at depinning (see also |13). 
We also note that in simulations of a rapidly moving 2D 
vortex lattice |5Sj the same nucleation mechanism as in 
Fig. El but relative to the comoving frame, was identi- 
fied as source of velocity differences between chains. 

For incommensurate channels close to commensurabil- 
ity, the TSF's which are caused by the mismatch are 
pinned by the disorder. For the weak disorder strengths 
considered here the random stress from the CE's is not 
sufficient to break up the TSF's. Consequently, at zero 
drive the array in the channel consists of a weakly disor- 
dered superlattice of TSF's. The behavior of the trans- 
port curves is shown in Fig. EJf° r w/bo — 3.03 and 3.06. 
It reveals features very similar to the curves of a single 
chain (Fig. I12|l . For / s < / < fi/n, a low mobility regime 
in the v-f curves develops due to the motion of the TSF's. 
Here f s corresponds either to the rms pinning force of in- 
dividual TSF's or to the collective pinning force for larger 
stack density. For forces / > ji/n the curves approach 
linear flow behavior again. 
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FIG. 19: (a): Time evolution of longitudinal displacements 
Ui.j of individual rows j — 1,2,3 at depinning for w/bo = 3 
and A = 0.05. (b) Square lattice representation of the nu- 
cleated stacks of discommensurations. Small arrows indicate 
the Burgers vector of the dislocations terminating each stack. 
The large arrows indicate their propagation direction, (c) 
vortex trajectories during nucleation of the vacancy stack be- 
tween x = 40ao and x = 65ao (up to the time corresponding 
to the filled symbols in (a)). 




FIG. 20: Critical force versus channel width for a disorder 
strength A = 0.05 using a velocity criterion v ~ O.Ol/i/7. 
The dashed line represents the continuum result f s = y,bo/w. 



B. Behavior around 'half filling' 

The dependence of f s on the channel width for A = 
0.05 is shown in Fig. |20l for two disorder realizations. 
For larger channel widths, smaller channel lengths were 
used with L > 1000 /(w/bo). The data around matching 
(w/bo ~ n) reflect the nucleation and pinning of TSF's as 
discussed above: f s at commensurability is reduced com- 
pared to the pure value (i/n and the commensurability 
peak is considerably broadened, particularly for larger n, 
due to the pinning of TSF's. The apparent discontinuity 
in the peak may be an artifact of the finite channel length 
(see the discussion in Sec lVBI on the distribution of nucle- 
ation sites). A new, and robust feature, however, is that 
around 'half filling' distinct maxima in f s appear. The 
origin of these maxima is illustrated by the static struc- 
ture near half filling. An example is given for w/bo = 3.48 
in for Fig. [2] The triangulation shows that, in addition 
to aligned dislocations with b || x, also misaligned dislo- 
cations appear with Burgers vector at an angle of about 
±60° with the CE's. These misaligned defects are lo- 
cally stabilized by the disorder in the CE's and thus also 
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FIG. 21: Triangulation of the static ground state structure for 
w/bo = 3.48 and A = 0.05. The arrows indicate the Burgers 
vectors of the dislocations. Shown is a characteristic segment 
of the total channel length L — 500ao. 
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FIG. 22: Vortex trajectories for the channel segment shown 
in Fig. |3T](«;/6o = 3.48, A = 0.05) during motion over ~ 4a 
at (a) / = 0.0038, (b) / = 0.0101, (c) / = 0.0209 and (d) 
/ = 0.0419. 

pinned by the disorder which leads to the increase in the 
threshold force. In addition, the projection of the driv- 
ing force along the glide direction is always smaller for 
misaligned dislocations than for aligned dislocations. It 
is seen that the misaligned dislocations separate regions 
with n rows from regions with n ± 1 rows. Either of these 
regions may thus be considered as a longitudinal stack- 
ing fault (LSF). In Sec II VI we mentioned that in absence 
of disorder LSF's are metastable, the structure with a 
single integer number of rows and a regular distribution 
of TSF's is energetically somewhat more favorable. In 
presence of weak disorder, however, the LSF are stable 
and have the lowest energy. 

As for the dynamics at / > f s , the dislocation struc- 
ture and flow pattern are generally different from the 
static pattern of LSF's. We illustrate this in Fig. |23 
where the vortex trajectories at various (increasing) driv- 
ing forces are shown. For small drive (Fig. I22f a). 
/ ~ 0.25/i / (w/bo)), a region of plastic motion within the 
channel is seen at about the location where the static pat- 
tern shows a 4 row structure. Vortex transport through 
these fault zones occurs by repeated nucleation and anni- 
hilation of misaligned defects. At a fixed driving force an 
LSF remains at the same position, although its bound- 
aries fluctuate over a distance of at most 2 — 3ao. This 
contrasts the situation in absence of disorder where LSF's 
can move along the channel via a 'climb' like process (see 

H). 

For different driving forces, the location and amount 
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of either norni 1-row regions or fault zones is differ- 
ent, as shown in Figs. (b),(c) and (d). In this particular 
segment the n = 4 region expands on increasing drive 
but at other locations the reverse can occur. Moreover, 
after cycling the force, a different structure can occur 
at the same drive, i.e. no unique structure exists at a 
given force. This may also lead to small hysteresis in the 
v-f curves. Overall, we see that the transverse degrees 
of freedom in the channel, in combination with disor- 
der, give rise to an important new mechanism of yield 
strength enhancement. In the following section we an- 
alyze these structures in more detail in the context of 
strong disorder. 



VII. WIDE CHANNELS WITH STRONG 
DISORDER 

The behavior of f s versus w/bo in Fig. [201 still exhibits 
considerable discrepancies with the experimental data in 
Fig. [21 Clearly, the CE disorder underlying these exper- 
imental data (see |59j) differs from the type of weak CE 
disorder considered so far. Motivated by recent imaging 
experiments |34j . which showed glassy vortex configura- 
tions in the NbN edge material, we now consider the case 
of strong CE disorder. As will be shown, in this regime 
the effect of transverse degrees of freedom and the pres- 
ence of misaligned defects provide the main mechanism 
for the critical current oscillations. 

The simulations we discuss here in detail were per- 
formed at a large disorder strength A = 0.2. The only 
remaining order in the CE arrays is their preferred ori- 
entation with the principle lattice vector along the CE's. 
The system sizes were typically such that wL > 1500a &o- 
We also allowed for quenched defects between the two 
CE's (uncorrelated longitudinal strains in the upper and 
lower CE arrays), but at these large disorder strengths 
this is not essential. 



A. Static structures, yield strength and depinning 

In Fig. |231 we show triangulations of the static vortex 
configuration for channels of width w/bo = 3, w/bo = 3.5 
and w/bo = 4. As seen in Figs, (a) and (c), even in the 
matching case dislocations are present due to the large 
random stresses from the CE. While some of these may 
originate from quenched 'phase slips' between upper and 
lower CE's, we checked that, also without these 'phase 
slips', the matching structure at this disorder strength 
always exhibits twisted bonds (defined by two adjacent 
dislocations of opposite charge) at the CE's as well as 
oppositely 'charged' TSF's terminated by defect pairs at 
the CE. While most defects are thus situated along the 
CE with b || x, occasionally they can be located inside 
the channel and have misoriented Burgers vector. Turn- 
ing to the mismatch case in (b), it is seen that a region 
of 4 rows (left) coexists with a 3 row region (right). In 




FIG. 23: Delauney triangulation of the ground state struc- 
ture at A = 0.2 for (a) w/bo = 3, (b) w/bo = 3.5 and (c) 
w/bo = 4. The arrows in (a) indicate the Burgers vectors of 
the dislocations. The aligned dislocations in (a) and (c) in- 
clude a disorder induced static dislocation pair (twisted bond) 
at the lower CE, in (b) numerous misaligned dislocations are 
present. 




FIG. 24: The 'correlation' length (L//) (o) and the fraction 
Li°i /L (•) of regions with n — 3 rows versus w/bo. 



between these regions, there must be dislocations with 
misaligned Burgers vectors. In addition, numerous other 
misaligned defects are visible, rendering local destruction 
of the chain alignment with the CE's (although less fre- 
quent, the latter can also occur in the matching state, 
see Fig. (c)). 

To further characterize the disorder in the / = struc- 
tures, we analyzed, for the regime 2.5 < w/bo < 3.5, 
the average length of domains without misaligned dis- 
locations, (£//), as well as the total fraction, Lj'f/L, 
of regions with n = 3 rows. The results are shown in 
Fig. 1241 As observed, both quantities are maximum at 
w/bo = 3 and decay considerably away from matching: 
for \w/b — 3| > 0.3 the average length of 'correlated' 
n = 3 regions in the static structure is no more than 
10ao and they make up less than ~ 50% of the channel. 
The remaining fraction 1 — (L 1 ?} /V) contains misaligned 
dislocations and small regions with 2 (w/bo — 2.5) or 4 
(for w/bo — 3.5) aligned chains. 

In Fig. [23 we show the behavior of the threshold 
force for A = 0.2. The modulation of f s with chan- 
nel width is still present but it has changed consider- 
ably compared to the weak disorder case. The sharp 
maxima at integer w/bo have vanished, very similar to 
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FIG. 25: Computed threshold force versus channel width for 
strong disorder, A = 0.2. The data were obtained by taking 
the friction force / — *yv at a velocity criterion v ~ 0.025(^/7) 
and subsequent averaging and smoothing over 5 disorder re- 
alizations. Drawn line: continuum result, f a = n/(w/bo). 



the case of the ID chain at strong disorder, see Fig. 
IPol Instead, we now observe smooth oscillations, with 
maxima in f s for w/bo ~ n + 0.65 and minima for 
w/bo ~ n + 0.15. The maxima, although occurring 
slightly above 'half filling', are of similar nature as the 
local maxima at w/bo — n±l/2 for weak disorder: they 
are related to numerous misaligned defects present in the 
structure around mismatch. They enhance the flow stress 
compared to that of the structures around matching with 
predominantly aligned defects. 

The differences in threshold force are also reflected by 
the vortex trajectories at small velocity. In Fig. 1261 we 
show these trajectories for channel widths w/bo = 3.1, 
3.6 and 4.1, close to the extrema in f s . The first thing 
to notice is that the trajectories at mismatch (Fig l26f b % ). 
w/bo = 3.6) are densely interconnecting, on a scale ~ ao, 
i.e. the motion is fully plastic and creation and annihi- 
lation of misaligned defects occurs over nearly the full 
channel length. For the 'near matching' cases in (a) 
and (c), the motion occurs mainly in the form of inte- 
ger chains. However, for the small velocity considered 
here, the dynamics still exhibits a considerable amount of 
plastic motion. This partly occurs due to vortices which 
remain stuck at the CE's (also visible in (b)) and partly 
due to narrow interconnecting regions. In fact, compar- 
ing the average length of the regions without inter row 
switching in (a) with the data for the static structure in 
Fig. 121 it is seen that near matching the structure at 
small v is more disordered than the corresponding static 
structure, reflecting nucleation of misaligned defects in 
regions which were free of such defects at / = 0. 



B. Analysis of dynamic properties 

We now consider in more detail the properties of the 
moving structures at / > f s . We first show in Fig. [TTT al 
two characteristic v-f curves associated with a minimum 
(w/bo = 3.1) and a maximum (w/bo = 3-6) in flow stress. 
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FIG. 26: Vortex trajectories at small velocity v ~ 0.001 for: 
(a) w/bo = 3.1 (/ = 0.0053 ~ O.lfi), (b) w/b = 3.6 (/ = 
0.0074 ~ 0.14/*) and (c) w/b = 4.1 ( / = 0.0044 ~ 0.08/t), 
all during motion over ~ 3ao- 
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FIG. 27: (a) v-f curves around a minimum in flow stress 
w/bo = 3.1 (o), and around a maximum w/bo = 3.6 (•). 
Inset: W x (xm) and W v (xm) for w/bo = 3.6 and / = 0.021. 
The dashed line emphasizes the linear behavior of W y in the 
regime 3ao < xm < IO00. (b) The density of switching points 
pi® = AWy/(AxM(io) versus force for w/bo = 3.1 (o) and 
w/bo = 3.6 (•). Inset: W x (xm) and W v (xm) for w/bo = 3.1 
and / = 0.041. The arrow indicates the value W VtC referred 
to in the text. 



It is observed that the enhanced threshold in the latter 
case also translates in a larger dynamic friction, / — jv, 
of the driven structure. In addition, the latter curve 
exhibits a small positive curvature in the small veloc- 
ity regime v < 0.01. This nonlinearity is related to the 
strong plastic nature of the motion in this regime. 

A convenient way to characterize the amount of plas- 
ticity is to calculate the mean square displacement of 
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vortices form their center of mass (M) positions 60] : 

W a {t) = + *o) - *i(t )] 2 /N ch , (39) 

i 

where a — x—xm, U~Vm denotes longitudinal and trans- 
verse displacements, respectively. As shown by Kolton et 
al for a 'bulk' 2D vortex lattice [§3], W a (t) can be char- 
acterized by W a = Rat^" . For example, when £ y = 1 
we have normal transverse diffusion (caused by 'random' 
switching of vortices between chains) with R y = D y the 
diffusion coefficient. However, in the channels W y will be- 
come bounded at long times (large Xm) due to the finite 
channel width. In the inset (i) to Fig. I27f a) we have illus- 
trated this behavior for w/bo = 3.6 and / = 0.021. For 
xm ^ 15ao (the point xm — was chosen in the steady 
state after transients had disappeared), W y increases lin- 
early as in usual diffusion, but for larger displacements 
(times) W y levels off and eventually saturates. In addi- 
tion, even when all vortices remain in their chain (no 
transverse diffusion), W y initially increases to a value 
W y>c due to finite chain 'roughness'. Such behavior is 
observed for xm ^ 2a in the inset (ii) to Fig. !27T b1 
where W y is shown for a more coherent flow situation 
at w/bo = 3.1 and / = 0.041. In practice, we found 
that W y<c is always reached for xm < 3ao, while for the 
long time (large distance) behavior significant levelling 
of W y occurs only when xm ^ IOoq. The appropriate 
regime we use to characterize real diffusion is therefore 
given by 3ao < xm < 10ao (i.e. Axm — 7ao). Further, 
D y itself does not directly reflect the density of chain- 
switching points along the channel. Defining the 'ID' 
density of such switching points as pi® ~ N sw /L, the 
rate of switching increases linearly both with pl° and 
with the average velocity: ijj ~ p\®v. Hence, the diffu- 
sion constant is given by: D y ~ a^/t sw — pl^va^. Being 
interested in pi®, we therefore divide out the intrinsic 
velocity dependence of D y and calculate W y (t) / (vta^) = 
&W y l(&x M al) = p l °. 

The results of p\® versus / are shown in Fig. I27f b) 
for the two cases w/bo — 3.1 and w/bo = 3.6. Around 
the minimum in flow stress (o), the overall value of p\® 
is considerably smaller than around the maximum (•), 
similar to what is seen in the trajectories in Fig. 1261 For 
both cases, is clearly reduced on increasing the force. 
This reflects both a decrease of the number of fault zones 
in the moving structure as well as a suppression of switch 
events (called 'transverse freezing' in 60] ) within regions 
already organized in n moving rows. Around matching, 
pi® smoothly vanishes at / ~ p ~ 0.05, indicating com- 
plete dynamic ordering into an n = 3 row structure with- 
out transverse wandering. For w/bo = 3.6 an ordering 
transition is also observed but it occurs at much larger 
drive (/ ~ 3/i, not shown) and the array orders into an 
n = 4 row configuration, with a reduced spacing b < bo 
between the chains and average vortex spacing a > ao 
within the chains. At the end of this section we will 
show how in this large drive regime the number of rows 
changes with w/bo- 



(a) 




FIG. 28: (a) Flow trajectories at large drive (a) during motion 
over ~ 3ao for w/bo = 3.1 and / ~ 2p. (b) Delauney triangu- 
lation of one snapshot of (a), (c) Trajectories for w/bo = 3.55 
and / ~ 2p. (d) Triangulation of a snapshot of (c). 



The insets to Fig. |57|also show the longitudinal mean 
square displacements. For strongly plastic flow at mis- 
match (inset (i)), W x is large and £ x is close to 2, as 
expected when some vortices remain stuck at the CE's. 
For the more coherent situation in (ii), where transverse 
switching has nearly ceased, W x is smaller and £ x > 1. 
We always find an exponent 1 < £ x < 2, similar to the 
results for 2D VL's in 60]. Interestingly, even without 
transverse wandering (p\® = 0), W x increases indefi- 
nitely (with £ > 1), indicating that the moving integer 
chain structure still exhibits slip events and (local) ve- 
locity differences between the chains. 

We illustrate some more aspects of the structures at 
large drive in Fig. GEJ where the vortex trajectories 
and triangulations of a single snapshot are displayed for 
w/bo = 3.1 and w/bo — 3.55, both at / ~ 2p. For 
w/bo = 3.1, we do not see any transverse wandering in 
Fig. EHta). Figure E^b) shows that in this case the 
dynamic structure exhibits only dislocations with b \\ x, 
mainly located at the CE but also occasionally between 
chains inside the channel. The latter dislocations are pos- 
sibly dynamically nucleated. They are non-stationary 
in the comoving frame and lead to slip events and the 
growth of W x as was discussed above. Turning to the mis- 
match case (Figs. I28f c) and (d)), it is seen that the dy- 
namic structure consists of n — 3 and n — 4 row regions 
coexisting in the channel. At the driving force considered 
here, p),® has a finite but small value a pl® — 0.007, 
which is solely due to switching of vortices in the fault 
zones separating the 3 and 4-row regions. Within these 
regions transverse wandering is absent. At yet larger 
forces the minority 3-row regions vanish and complete 
ordering into 4 rows occurs, as for the case w/bo = 3.6. 
For a given driving force, the n-row regions again remain 
quasi-static during motion. The triangulation in (d) ex- 
hibits the expected misaligned dislocations at the fault 
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zone but in general also aligned dislocations are present 
between the chains within an n-row region (not shown in 
the figure). 

It is also interesting to compare the velocities in the 
two coexisting regions with n and n ± 1 rows. Denot- 
ing the vortex velocity in an n-row region by v n and 
the longitudinal vortex spacing there by a n , flux con- 
servation implies that nv n /a n = n'v n >/a n i. We checked 
that, in both regions, the average flux density \/{a n b n ), 
with b n — w/n the row spacing, was equal to l/(ao&o) 
within ~ 4%. Therefore a n / — (n'/n)a n , and conse- 
quently the average vortex velocities are equal, v n = v n i . 
However, the local washboard frequency, v n = v n /a n , is 
different in both regions. Indeed, the spectrum of the ve- 
locity fluctuations in channels with dynamic coexistence 
of n and nil rows showed two shallow fundamental 
peaks at frequencies v n /v n l — n'/n. We however note 
that, both around matching (where a single peak occurs 
at v ~ v/clq) and around mismatch, the amplitude of 
the washboard peak(s) decays on increasing the chan- 
nel length L. In addition, for large velocities the mixed 
n/n i 1 structures ultimately anneal into a single n or 
nil domain, causing the collapse of one of the peaks. 

The simulations also allow to explicitly show the in- 
fluence of the transverse degrees of freedom on the mod- 
ulations of the dynamic friction force (and, ultimately, 
the critical current). Generalizing the expression for the 
friction force Ea. (|20|l in Sec. IIII Bl including transverse 
fluctuations leads to: 

/ =( 7 /v)[{(d t u x f) itt + ((d t u y f) i , t ] 
= 7v + ff ric + fj ric 

= 1V + f fric , (40) 

where /?_,•_ and fU- denote the contribution to the total 

J i tC "J VIC 

friction due to longitudinal and transverse fluctuations, 
respectively. Figure [5^ a) displays the numerical results 
for these quantities, normalized by the total friction force 
/ — l v as obtained from the v — f curves. The sum of 
the two data is 1 as it should be, confirming the correct 
numerical evaluation of these quantities. The total fric- 
tion is also shown for clarity (Fig.(b)) and is essentially 
the same as the data shown in Fig. [23] Clearly, the 
relative contribution of longitudinal fluctuations to ffric 
decreases on approaching a mismatch situation, while 
ffric/ ffric increases accordingly. These qualitative fea- 
tures remain present also at larger velocities, where per- 
meation modes between chains are being suppressed. 

Finally, we consider the continuous modulation of 
structural (dis)order in the moving arrays when vary- 
ing the channel width. As was shown in Fig. I27f b) , for a 
fixed velocity, the density of switching points p\® is max- 
imum around mismatch while for fixed channel width it 
decreases with velocity. We can then define an order- 
ing velocity v c operationally as the point at which p\® is 
reduced below a certain threshold. In Fig. I29f c~) the be- 
havior of v c for channel widths w/bo > 2 is shown for two 
criteria. The lower curve (v c ,i), with aop\® w 0.01, cor- 
responds to a situation with ~ 70% of the channel length 
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FIG. 29: (a) The contributions ff ri J f fric (•) and ff ri Jff ric 
(o) to the dynamic friction versus w/bo for a velocity v = 
O.l/i/7. (b) Total friction / — at v = O.I/1/7. (c) Ordering 
velocity v c ,i (•) determined from a criterion aopl® — 0.01, 
and v c ,2 (o) using a criterion aopl® — 0.002 (error bars were 
estimated from different disorder realizations), (d) Number 
of moving rows n at large drive (v > v c ,2) versus w/bo- 

'transversely frozen' into integer chain regions of length 
> 10a . The upper curve (v c ,2), with a pl® « 0.002, 
corresponds to nearly fully annealed arrays. As observed, 
f c ,i increases smoothly by about an order of magnitude 
between a matching and mismatching situation, while 
v c _ 2 increases by a factor > 5 (except for w/b > 5). Re- 
gardless of the criterion, the amplitude of the modulation 
of v c is considerably larger than the amplitude of the f s 
modulation, which we will further discuss shortly. 

For v > t» Cj 2 the arrays all completely anneal into a sin- 
gle n-row structure without permeation modes. Figure 
l2*9T d) displays the number of rows n of these structures 
versus w/bo. The switching from n to n + 1 rows is seen 
to occur at half integer channel widths for w/bo > 3 but 
the transitions 1 — ► 2 and 2 — > 3 take place below these 
points, at w/b ~ 1.35 and w/b ~ 2.4, respectively. 
Around the transition points, there are regions in which 
n and nil rows coexist at smaller velocity v < v Cj 2- The 
widths of these regions were found to be Aw/bo — ±0.05. 



VIII. DISCUSSION 

The results for f s versus w/bo at strong disorder in the 
previous section show strong resemblance to the measure- 
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FIG. 30: Comparison between the experimentally measured 
yield strength of an artificial flow channel and the numerical 
results at strong disorder (rilll , represented through the param- 
eter A (see the text) normalized by the value A = (7r\/3) -1 
for the ideal lattice. 



merits of the critical current versus field B in the exper- 
imental channel system, see Fig. |5J For a more detailed 
comparison, we combine in Fig. 1301 both data, repre- 
sented in terms of the parameter A (see Sec|I]): for the 
experiment, A exp is determined from A exp = F s w / (2cqq) 
and is plotted versus B 1 / 2 , for the simulations, A — 
(f s /fi)(w/bo)A° As observed, the shapes of the 

oscillations are in very reasonable agreement, although 
the overall value of A from the numerics is larger than 
the experimental value [6^ | . The important conclusion 
from the simulations here is thus that the maxima in 
the measured critical current do not correspond to tra- 
ditional commensurability peaks, but are caused by en- 
hanced plastic motion and transverse deviations in mis- 
matching channels with strong edge disorder. Previously, 
we indeed obtained experimental evidence for this mecha- 
nism via mode-locking (ML) experiments 31] in which an 
rf-drive is superimposed on the dc-drive. Interference be- 
tween the former and collective modes of deformation in 
the moving array leads to plateaus in the current-voltage 
curves. The plateau (ML) voltage then directly yields the 
number of moving rows n and transitions from n — > n± 1 
rows were observed to coincide with maxima in critical 
current 1631. 



More detailed ML-experiments provided a wealth of 
additional information on the dynamics of the arrays 
|32|. In particular, a minimum velocity was required 
to observe the ML-phcnomenon, which we identified as 
ordering velocity. This ordering velocity exhibited a 
strong upturn away from matching, similar to the be- 
havior of v c in the simulations. As proposed in |32j |. 
v c can be estimated from a modified version of the dy- 
namic ordering theory of Koshelev and Vinokur (KV) 
[37j . Firstly, instead of the 2D random potential in (37) , 
for the channels v c is related to the short range corre- 
lated random stress from the CE, with r.m.s. amplitude 
~ e ce C66 and e ce the random strain. Secondly, when ther- 
mal fluctuations can be ignored compared to typical de- 
fect energies (see below), v c is inversely proportional to 



the energy ksT p for creation of small defect pairs (see 
also H3|)i v c /a ~ y / 3/32Tr(a /a) 2 (e ce c 66 ) 2 aobo/(-/kBf p ), 
with kgTp ~ A p CQQa 2 l /2ir (per unit vortex length). The 
typical pairs referred to are those with misaligned Burg- 
ers vector, which are responsible for breakup of the chain 
structure. Hence, the increase of v c away from matching 
implies a decrease in the pair formation energy, which 
is accounted for by including the parameter A p < 1 in 
ksTp, while setting A p = 1 at matching. Experimentally, 
A p decreased to ~ 0.1 at mismatch. The behavior of v c 
and A p in the simulations is analyzed using the dimen- 
sionless form of the above formula for v r : 



(7« c /m) - 2.5e 2 ce (a /a)/A p 



(41) 



The data for u Cjl around matching (where A p = 1) yield 
as measure for the random strain e ce ~ 0.19 and, near 
mismatch, a reduction of the defect pair creation energy 
by a factor A p ~ 0.1. The latter is in very reasonable 
agreement with the experiments, and the reduction of 
ksTp near mismatch also qualitatively agrees with the 
large number of disorder induced fault zones observed in 
the static structures near mismatch (Fier 123(1 . However, 
the result for e ce is considerably larger than the value 
£ce = 0.025 found experimentally. This is also manifest 
in the fact that in the experiments the pinning frequency, 
defined as / s /(7»), always exceeds the ordering frequency 
v c /a, while in the simulations f s < jv c . 

Within the modified KV theory, the pinning frequ ency 
and the ordering frequency are directly related via |32l |: 
v c /a = T(f s / ia ) 2 with t = (e 2 ce /2A 2 )(wB) 2 /(ApC 66Pf ). 
Experimentally, r was found to be independent of the 
matching condition. The decrease in the defect pair cre- 
ation energy (oc A p ) then relates to the increase in the 
pinning frequency (yield strength) away from matching 
(oc A) as A p tx 1/A 2 . The numerical results for A p versus 
A behave similarly and can be fitted by A p tx 1/A S with 
<; ~ 2 — 3, but data collected over the full range of w/b 
show to much scatter to make a more detailed compari- 
son. Nevertheless, the KV model qualitatively accounts 
for the enhanced amplitude of the v c modulation com- 
pared to that of the yield strength. 

As for the discrepancy between the experimental and 
numerical values for e ce or v c , one should keep in mind 
that experimentally the ordering velocity is determined 
from the onset of an n-row ML plateau, i.e. it corre- 
sponds to the velocity at which coherent n-row regions 
first appear, while incoherent regions may still exist in 
other parts of the channel. In the simulations this may 
thus correspond to v c determined using a larger criterion 
for the density of switching points p\®. In addition, the 
superimposed rf-drive in the experiments may assist re- 
ordering of the structure, also leading to smaller values 
of v c . For future studies it is interesting to directly com- 
pare numerical simulations of channels with mixed rf-dc 
drive with the experiments and test which criterion best 
represents the reordering phenomenon. 

Additionally, the experimental results indicated that 
in the large drive regime (where the ML amplitude satu- 
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rates at a constant value) , the coherently moving fraction 
of vortices does not exceed ~ 40% (at matching) , while 
it was reduced on approaching mismatch. This feature 
appears at odds with the simulations where eventually 
the arrays all order into a completely transversely frozen 
n-row domain at large drive, regardless of the matching 
condition. At present we do not have a good explanation 
for this discrepancy. 

Finally we comment on the T = approximation made 
throughout this study. When thermal fluctuations are 
important, not only does one expect the dynamic or- 
dering behavior to be changed (see [33, EU), also the 
sharp threshold behavior in the v-f curves will be smeared 
and activated flow may occur. The relevant energy scale 
for both phenomena is again the energy for formation of 
small defect pairs ksT p . In the low magnetic field exper- 
iments in [23| this energy was ~ 2 orders of magnitude 
larger than UbT and for comparison with these results 
the T = approximation is justified. To compare with 
the experiments near the melting field in |65j |. it would 
be required to include thermal fluctuations in the simu- 
lations. 



IX. SUMMARY 

We have presented a detailed study of the properties 
of vortices confined to narrow flow channels with pinned 
vortices in the channel edges. In the experimental sys- 
tem which motivated this work [2^ . l3ll l32| , the thresh- 
old force (yield strength) shows pronounced commensu- 
rability oscillations when the natural vortex row spac- 
ing is varied through integer fractions of the channel 
width. The analysis and simulations presented in this 
paper show that in a mesoscopic channel system the de- 
pendence of threshold on commensurability as well as 
the dynamics of vortices in the channels drastically vary 
with the amount of disorder in the confining arrays. At 
zero or weak disorder the system behaves similar to ID 
LJJ-systems and defects at the CE's reduce the yield 
strength. At large disorder the behavior involves transi- 
tions from quasi-lD to 2D structures, where an increase 
in the amount of plastic deformations enhances the yield 
strength similar to the situation in the classical peak ef- 
fect in superconductors. 

We first presented a generalized sine-Gordon descrip- 
tion for a ID vortex chain in an ideally ordered channel. 
In this case (or for channels with multiple chains near 
commensurability) the threshold force has sharp peaks at 
commensurate widths, whereas it is essentially vanishing 
at incommensurability due to easy glide of 'aligned' de- 
fects, i.e. defects with Burgers vector along the CE. The 
model was then extended to study the effects of weak 
disorder in the confining arrays. Simulations and ana- 
lytical results showed that this disorder causes the sharp 
maxima in the threshold force at matching to be lowered 
and broadened due to nucleation (at matching) and pin- 
ning(at mismatch) of edge defects. Apart from these de- 



fects, the arrays respond elastically in this regime, both 
near threshold and at large drive. We studied numeri- 
cally the relevant edge defect dynamics and, using the 
sine-Gordon model, we analyzed the crossover to strong 
defect pinning on increasing the disorder strength. 

For large disorder in the CE's, matching between the 
longitudinal vortex spacings in and outside the chan- 
nel becomes irrelevant and the peaks in threshold force 
around matching completely vanish with a 'saturated' 
value for f s of about 30% of the ideal lattice strength. 
Around mismatch however, the arrays become suscepti- 
ble to formation of defects with Burgers vector misori- 
ented with the channel direction. Such defects either lo- 
cally break up the integer chain structure or exist at the 
boundaries of n and n ± 1-row regions coexisting in the 
channel. At large disorder, they are strongly pinned and 
cause the threshold force to exceed that around match- 
ing. Approaching a matching condition, the density of 
misaligned defects is reduced and a smooth modulation 
of f s results, with minima near matching. The depinning 
transition always involves plastic deformations inside the 
channel, but the amount of plasticity drastically increases 
away from matching. Using the density of transverse 
switching points (obtained from the transverse diffusion 
in the moving structures) as dynamic 'order parameter', 
we study the evolution of the moving structures on chang- 
ing the channel width and the drive. The arrays reorder 
(partially) into transversely frozen n-row regions at a ve- 
locity v c which shows a similar modulation with commen- 
surability as the threshold force. Finally, we compared 
the modulations of f s and v c at strong disorder with the 
available experimental results and with the dynamic or- 
dering theory in 37| and find good qualitative agreement. 
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APPENDIX A: ORDERED CHANNEL FOR 
ARBITRARY FIELD 

In this appendix we calculate the edge potential and 
sine-Gordon parameters in a symmetric, ordered channel 
with w ~ bo for arbitrary field. The interaction between 
a vortex at r = (x,y) in the channel and the pinned 
vortices in the CE's at R„ !m is given by: 

V ce (r) = (2n)- 2 J dk]Tu(k) e ik -( r - R "-), (Al) 

where V(k) is the Fourier transform of the vortex- vortex 
interaction. To obtain an expression valid over larger 
range of fields than the (low field) London-regime, we use 
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a generalization of the London potential Eq.© as pro- 
posed by Brandt [3(| • This generalization accounts for 
the reduction of superfluid density with field and an ad- 
ditional attractive interaction due to overlapping vortex 
cores. The Fourier transform of this interaction reads: 



V(k) = 2irU Q (l - b) 



1 



1 



|k| 2 + A' 



,(A2) 



where b = B/B c2 , A' = A/(l - b) 1 / 2 and £' = g (2 _g )1/2 
the effective coherence length with C ~ 1. It is conve- 
nient to split Eo. (|Al|) for the total potential in terms of 
the contribution V m of row m (see Fig. ©■ Integrating 
over k y and using Poisson summation yields 

V m (r) = \Vi(y)\ ™slk (x - ma /2) (A3) 



formula of Brandt Eq. © . Additionally, the edge poten- 
tial Ea. HA5|) is harmonic for all fields. Hence the ideal 
flow stress of a commensurate, ordered channel, is char- 
acterized by Frenkels value A = 1/ 7r\/3 for all fields [6(| . 
In the low field limit b < 0.2 (and X/clq > l/7r) the above 
expressions reduce to Eq.©, © in Sec. lIlTAl 

Using the field-dependent vortex interaction Eq. l|A2(l . 
one can derive the parameters in the sine-Gordon de- 
scription of Sec. IIII Al generalized for higher field. The 
equation for the chain stiffness becomes: 



Kq = £/ 7r(l - b) 



A'/ao 



y/1 + A' 2 g 2 y/l + £'V 



(A8) 



The reduced stiffness is obtained from g(b) 
K (b)/(k n(b)al): 



where ko — 2n/ao and the prefactors are: 



\Vi(y)\ = k U (l-b) 



~ lk o,\>\y', n \ 



' lk o,c\y' m \ 



Ik 



0,A' 



lh 



OA' 



(A4) 



W ith lk ,y = ^/(lk ) 2 + (A')" 2 , lk , t 

V( lk o) 2 + (?)~ 2 and V'm = -y±_m[b + ((w-b )/2\m\ 



First of all, we neglect in Eq. (|A3(l the / = terms which 
represent uniform (x-independent) interaction. Sec- 
ondly, when A > ao, as practically encountered in 
films, fco.A' ~ fco- Then the \l\ > 1 terms can be 
neglected, resulting in a sinusoidal potential |^| = 1. 
Further, in the summation over m only the contribu- 
tions from the m = ±1 terms are significant. Next 
we employ the relation £ 2 /a 2 = 6\/3/47r and rewrite 



k ,e = (l/a )^/47r 2 + 8vr(l - b)/{y/iC 2 b). The resulting 
expression for the total edge potential is: 



V cefi {v) 



cosh(fc y) 



,(1- S (6))tt 



S (b) 

. , - H>Q 

x2J7q(1 — b)e a o cos/co^, 



■ cosh(k ys(b)) 
(A5) 



7(6) 



3tt 



V3b 



f(b)V 4tt(1-6) U V2 



A C 



(A9) 



Taking into account these refinements in Ea. (|10|l and 
Ea. (|17|) . the defect width in the nonlocal regime becomes: 



lf(b) = 6Tr 2 a Q /f(b). 



(A10) 



One can obtain the typical crossover field b n i (or typical 
A/ao, ra i) at which nonlocal behavior sets in for a chain 
in an ordered channel by equating Eo. (|A10|l to the s- 
G value for the kink width 2nao \fg{5) ■ Approximating 
f(b) ~ (1 — 6), one finds 1 — b = 3tt/ ^g(b), which has 
the approximate solution 



Vnl — 77 



1- Jl- (48%/37r 3 e 2 /A 2 ) 



(All) 



Hence, the nonlocal regime is absent for a channel in 
a material with A/£ < 50 (and thickness d > A). For 
V£ ~ 60, nonlocal behavior occurs for b > b n t with 
bni ^ 0.2. This is to be compared with the estimate 
oo.ni < A/37T resulting from a simple London interaction, 
see Sec. IhTaI 



where e(b) = y/l + (2- 2b)/(nVSC 2 b). 

For a channel with w ~ bo, the maximum /i(b) of the 
sinusoidal pinning force —d x V ce o at y = is then given 
by: 



Kb) 



,C/o(l-6)/(6) 
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where we used e _7rv ^ ~ 1/247T 2 and 

/ Tv~/3{2-2b) 

„ P rv — c^i — 

f(b) = 1 - 24^ 3 



7rg-26) 2 



(A6) 



(A7) 



It can be checked that the associated shear modulus 
c 66 = 7rv / 3/i(^)/2ao is very similar to the interpolation 



APPENDIX B: SOLUTION TO THE DYNAMIC 
SINE-GORDON EQUATION 

For a displacement field of the form Eq. 124|) , expressed 
in modes with wave vector mq = 2irmcd and amplitude 
h m , Ea. H23|) for the v-f curve attains the form: 



f = yv 



2k 



M 



rn 



(Bl) 



The amplitudes h m are obtained by inserting Eq. I|24|) into 
the equation of motion l|13l) with the wavelength depen- 
dent elasticity n(q) from Eo. (|12|l . Since h can become of 
the order of a lattice spacing ao, one expands the 'sin' 
term in l|13fl up to second order in h. Furthermore, we 
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assume that h m decays rapidly upon increasing to and we 
keep only the three lowest order contributions h m with 
to < 3. Collecting terms of equal wave number, one ob- 
tains the following set of approximate equations: 



{ijvk Q + K hq )hi = -i + itfdhil 2 + \h 2 \ 2 ) 

+k h 2 - i{kl/2)h\ (B2) 



we rewrite Eo. (|B9l 



k v* = K eff (c d )/j. 



(Bll) 



Here r y/K e ff(cd) can be interpreted as the effective relax- 
ation time for the non linear deformations in the chain, 
which is expressed through the relaxation times ~//K m ^ q 
of linear modes (phonons) by Ea. (|B10|) . The velocity 
dependence of ^2(m\h m \) 2 may then be written as: 



— (2ijvk + K 2 , q )h 2 — kohi + ik Q {h\h 2 — hih 2 ) (B3) 



/' 



{3ijvk + K 3iq )h 3 = k h 2 + i{kl/2)h\, (B4) 



where K„ hq — m 2 q 2 K(mq) and h* denotes the complex 
conjugate of h. At small v, the real components of h 
vanish and the amplitudes describing the shape of the 
quasi-static deformations are given by: 



+ ^ (B5) 



3k 2 



-ko\hi 



2k 2 \h l \-{2K 2 Jn) 



Mh 2 \-kl\hx\ 2 /2) 



2K, 



3,9 



(B6) 



(B7) 



For arbitrary v, the /i m 's in Ea. (|B2IB3IB4(l may be deter- 
mined by a mathematical program. However, the most 
important effect of the coupling is that above a char- 
acteristic velocity v* (see below), \h 2 \ and \ha\ decrease 
rapidly as \h 2 \ ~ j/f and \h 3 \ ~ \hi\ 2 /3v, so that only 
the first Fourier mode survives. This mode is described 
by (v > v*): 



hi{v) = 



(2-yv/fi) - y/(2 7 ^//i) 2 

fc 



2jvko 



(B8) 



Since Ea. (|B8|l can be obtained by neglecting the elas- 
tic force terms K n ,qj it describes the deviation from the 
average velocity of a single particle in a periodic poten- 
tial. Further analysis of Ea. (|B2IB3IB4(l shows that the 
crossover velocity v* is determined by the amplitude of 
the quasi- static result in Ea.l |B5IB6IB7| l: 



(2 1 v*k /^) 2 J2( m \h m \) 2 ~l- 



(B9) 



Using the definition 



2n ) ~ 4 



^ (m\h„ 



, (BIO) 



E 



TO 2 I fir, 



A[KJ 



eff 



{jvk y 



(B12) 



This has the correct small v behavior (where higher 
modes play a role) and large v behavior (where h 2 w h 3 w 
0). First order perturbation using only h\ to first order 
in Ea. l|B2|) yields the same functional form but with K 2 ^ 
replaced by K\ , supporting the analytical interpolation 
made in obtaining Ea. (|B12|l . Finally, using Ea. lBl|) one 
arrives at Eq. l|25|l in section ITlI Bl 



APPENDIX C: ELASTIC SHEAR WAVES IN 
COMMENSURATE, ORDERED CHANNELS 

Starting from Eq. (|26|l we write u as the sum of a spa- 
tially uniform dc-component and a non-uniform, time- 
dependent component h(y,t). The result for h(y,t) — 
u(y,t) — vt is the following modified diffusion equation: 



1 d t h(y,t) = (f- 1 v)+cd 2 y h(y,t), 



(CI) 



where c = b^kQ^/2 = c^a^bQ. The boundary condition 
is set by: 

jd t h(-w/2, t) = -{fx/2) sm[iu Q t + h(-w/2, t)} 

+ (c/b 2 )Ah(t), (C2) 

and similarly for y = +w/2. Here the discrete char- 
acter of the array near the CE's is retained in Ah = 
h(—w/2 + bo, t) — h(—w/2, t). For the pinning term (last 
term on the r.h.s.) we now use the large velocity expan- 
sion for the restoring force from one CE: (fJ,/2) sin[fco(ui+ 
h)] ~ (n/2)sia[uot]. Then Eas.l(CT |l .l(C2 )l become similar 
to those for heat diffusion in a rod with at both ends 
heat sources that vary sinusoidally in time. In our case 
the frequency is the washboard frequency loq = k^v. By 
separation of variables one finds: 



h(y,t) = h e A etV (w)[f v (y) cos(oj t) + g v (y) sin(u> t)] 

(C3) 



(/-7«) 2 

-y , 



2c 

where f v (y) = -coa(y/l± tV )cosHa.(y/l± tV ) and g v (y) = 
sin(y/Zj_,t>) sinh(y//j_,„) with a velocity (frequency) de- 
pendent 'healing' length Zj_ -t , = y/2c/^ujQ = y/ (ppyv)bo. 
The factor A e , v (w) = 1/ y/f 2 {w/2) + g 2 {w/2) normal- 
izes the displacement at the CE to h e . The latter is ob- 
tained from the boundary condition Eq. I|C2|) . In the limit 
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t(Tfao z /lJb) 



ytaa 



lx,v < 2£>o- Since we have a discrete number of vortex 
rows, it is the two outer rows at y — ±w/2 which give the 
dominant contribution. Evaluating / — yv using Ea. ljC4|) 
then leads to Eq.® in Sec. HVI 



APPENDIX D: DISORDERED CHANNEL 
POTENTIAL AND PINNING OF DEFECTS 



FIG. 31: Velocity profile d t h(y,t) versus time from Ea. ljC3fl 
in a channel with w/bo — 9. (a) for 'yv/n ~ 0.4 and (b) for 
7^//i — 2. 



w/bo = n > 3 and yv/ p > 0.25, h e can be approximated 
by: 



2yv 



2kojv V 271; + p 



(C4) 



The solution (|C3|) describes a periodic velocity modu- 
lation dth of each chain, with a y-dependent amplitude 
\h'\ and phase shift. The latter reflects periodic lagging 
and advancing of the inner rows with respect to the outer 
ones. Both the decay of \h'\ away from the CE and the 
phase shift strongly increase with velocity through lj__ v . 
In Fig. we have illustrated this behavior for small 
velocity yv/p — 0.4 (a) and larger velocity yv/p = 2 
(b). Although the result in (a) is actually out of the 
regime of validity of the high velocity expansion, the 'in- 
phase' behavior at small v is a feature that qualitatively 
agrees with the simulations. For comparison we show 
in Fig. |23 numerical results for the oscillating velocity 
component for w/bo = 9 and small drive, / ~ 2f s ~ /x/5 
(yv/ p w 0.16). In addition to the small phase shift one 
observes that the modulation is highly non-sinusoidal, 
which is not captured by our approximate solution. 

Based on the above, one can also evaluate the dy- 
namic friction force f — yv. Using Ea. H20l) . 123|) . the 
expression for the elastic-continuum u(y,t) would be 
/ — yv — (y/vw) f dydt[d t h(y,t)] 2 . However, in the 
velocity regime where our solution applies, the length 



In this appendix we derive the disorder corrections to 
the channel potential and the effect on pinning of the 
chain. For the vortex interaction we use the London po- 
tential of Eq.©. Starting from Ea. (|30|l we first split the 
potential in contributions from row m with transverse 
coordinate y: 

V m {x) = (XT 1 J dkV(k, y) Pm (k)e lkx , (Dl) 

where k = k x and 



V(k,y) = U 7T- 



Vk 2 + A- 2 



(D2) 



In addition, p m (k) in Ea. (|Dl|) is related to the longitu- 
dinal displacement field d^x) via 0] 

p m (k) = a^ 1 / tto [1 - d x d m + 5 Pm (x - d m )] e~ ikx (D3) 



which contains the microscopic density modulation 
5p m (x - d m (x)) = J2i^o cos (l k o[x - d m (x)]) as well as 
the density variations cx d x d m due to long wavelength de- 
formations. Inserting Ea. (jD3jl in Ea. (|DlJI yields a total 
potential of the form V m (x) = Vi tm (x) + V p , m (x) where 

Vi, m (x) = -a^ 1 [ dx'V{x-x',y)d x ,d m {x'), (D4) 



and 



V p>m (x) = — / dx'V(x - x',y) cos[k (x' ~ dm)]- (D5) 
a 
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A constant offset arising from the 'DC'part of the den- 
sity has been omitted in Ea. (|D4|l and in Ea. (|D5|l only 
the I = ±1 components of the density are taken into 
account. The latter represents a quasi-periodic potential 
with wavelength ~ ao- It is only significant for y < 1.5&o, 
i.e. for m = ±1, as in the pure case. Therefore V p can 
be written as 



V p = —[jj, + Sp,(x)] cos[k (x - d)]/k , 



(D6) 



where d = (d\ + d_i)/2 and 

sn{x) =m dk[v(k + ,b Q ) 

~ nVSfidxd, 



V(k-,b Q )}ikd(k)e lkx 

(D7) 



where k± — k ± k and the second line is valid for 
k < 0.4fco- The local strains thus also provide varia- 
tions in the potential height. A similar conclusion holds 
when transverse shifts in the CE are included. In the 
simulations the mean square strain is ((V • d) 2 ) = A 2 /3 
and short range (~ ao/2) correlated along x with cor- 
relator S(s,0) = (d x d(x,y)d x d(x + s,y' = y)) x ~ 
(A 2 /3)e _ ( 2s / a °)~. Accordingly, the amplitude fluctua- 
tions of the periodic potential are characterized by: 



(6p,{x)6[i(x + .s)) 3 
k 2 



(/xAa /2) 2 e~ ( ^ ) (D8) 



The non-local contributions Eq. (|D4|1 add up to a total 
contribution Vi — X) m ^o m V ce . Vi thus originates 
from strains within a region ~ A around the channel and 
will be smooth on the scale ao (we assume A 3> ao). 
Transforming the sum into an integral and using Eq. (|D2H , 
the correlator Vi = (Vi(x)Vi(x+s)) can be written as |67| : 



a JCE 



dydy' / dxA Viyl (x)S{s - x, y' - J/)(D9) 



The term A y y(x) = {in)- 1 J dkV(k, y)V(k, y') cos(fcc) 
can be approximated by 



Ay,y' ( x ) 



\Jl-K\e 



lal+lsij 



-(f) 2 



(D10) 



In case the strains are uncorrelated between the rows 
S(x, y'-y)~ exp(-(2( a; 2 + (y' - y) 2 )/a 2 )) and Eq.^ 
can be approximated by 



T,( S ) ~C a A 2 U 2 (\/a ) 1+a e-^ 



(Dll) 



with a = 1 and C\ ~ 1. In case of uniform strains, S(x) 
is independent of y' — y. Then again the correlator is 
given by the above formula, but with a = 2 and C2 = 
(4/3 V 3 / 2 . 

To study the effect of the disorder on the pinning of 
vortices inside the channel, we write the total energy as 
a sum of an elastic and a pinning term: 



H{x, u) 



H e i + H p = a^ 1 J dx^(d x u) 2 + # P .(D12) 



The dispersion of the elastic constant has been neglected 
and u represents the displacements of vortices inside the 
channel. By writing the density of the chain as p c (x, u) = 
% 1 [l — d x u + '^2 l _ i0 cos(lko(x — u))], the pinning term H p 
111 can be expressed as: 



dx(Vp + V l )(dp(x,u)-d x u), (D13) 



where only the lowest Fourier components of p c are re- 
tained, i.e. Sp(x 1 u) = 2cos/co( :r — u )- This expres- 
sion can be simplified as follows. Since we consider the 
limit where d x u is nearly constant within ao, the cross- 
term of the quasi-periodic potential V p and d x u can be 
neglected compared to the Vi{x)d x u term. The prod- 
uct Vi6p(x, u) is also oscillatory and can be neglected as 
well |68| . The other remaining term V p (x)Sp(x,u) can 
be written as (fj,(x)/ko) cos[£;o(u — d)]. Shifting the ar- 
gument via k (u — d) = k u and writing u — ► u, an 
extra term K d x dd x u is generated in H e i in Ea. (|D12|) 
(and also u-independent terms which can be neglected). 
The total coupling to the strain d x u then consists of 
V s = Vi — Kod x d. The result for the total energy is 



H = H S g — I — 



dx 
ao 



5fx(x) 



cos(kou) — V s (x)d x u 



(D14) 



The potential V s has a correlator (V s (x)V s (x+s)) — T s (s) 
which is characterized by 

167rV 2U A 2 X ^ f 

r s( s ) = — o — r a(s) + F|(s) H V{s, 6 )(D15) 

3 a 

in which g — kq / (2n /iao) and the last term 
arises from cross-correlations. In Eq. ijD14fl . H§q = 



J dx[(K /2)(d x u) 2 — (n / k ) cos(k u)] is the original sine- 
Gordon energy functional of the pure model and the re- 
maining terms reflect the corrections due to disorder. Fi- 
nally, we denote the first correction, which is due to the 
amplitude fluctuations, as H a and the second, coupling 
to the strain, as H s . 

For weak disorder, we can now calculate the ef- 
fect of disorder on a defect in the chain by assum- 
ing that the shape of a defect at x, Ud(x' — x) — 
2ao arctan[exp(±27r(a;' — x)/ld)]/iT, is unchanged by dis- 
order 0| . The pinning energy of a defect due to the 
term H a is: 

E a {x) — (ao/co) 1 J dx' Sp,(x') cos(koUd(x' — x))(D16) 
The correlations of E a are given by |6^ |: 

{E a (x)E a (x + s)) = -j f d P A a (p)T a (s- P ), (D17) 



where 



Aa( P ) 



4(l d /2ir)dx 



cosh 2 (x) cosh 2 [x + p) 



ie ^ } , (D18) 
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with x = lux I Id- For Id 3> ao we can approximate T a (s) 
in Eq. ^nnt by r a (s) ~ (^A) 2 (0Fag/8)<5(s) leading to: 

(£ a (x)£ a (.T + S )) = (V^/8)/i 2 A 2 / d a e' ( ft )2 . (D19) 

The effect of the coupling to the strain is given by 
E s (x) = (ao) -1 J dx'V s (x')d x Ud{x' — x) which has the 
following correlator: 

(E s {x)E s (x + a)) = \ I d P A s {p)T s (s- P ). (D20) 
a o J 

Using d x u d = (2ao/^)cosh~ 1 (27rx/Z d ), A s is given by: 
(2al/irl d )dx 



Mp) 



4a 2 e~ { ^ ) 



cosh(i) cosh(ic + p) 



■ (D21) 



For > A the final result is: 



(£ s (z)£; s (x + s)) 



(t/pAA) 2 



[Q+ije^ , (D22) 



where the term C; ~ 2C a (X/ao) a in square brackets is 
due to the nonlocal contribution Vj and the factor 4 arises 
from the term ~ d x d in y s . Hence, for large A/ao the 
nonlocal term dominates in the coupling to the strain. 
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